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Abstract 

It  is  demonstrated  that  a  nonlinear  Schrbdinger  equation  with  certain  nonlineari¬ 
ties  allows  for  an  existence  of  multi-stable  singular  solitons  (i.e.,  singular  solitons  with 
the  same  carried  power  but  different  propagation  parameters).  In  nonlinear  optics, 
these  solitons  may  exist  either  in  the  form  of  short  bistable  pulses,  or  bistable  self¬ 
trapping  (both  two-  and  three-dimensional). 
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In  this  letter  we  demonstrate  that  for  a  certain  class  of  nonlinearities,  the  soliton 
solution  of  the  (generalized)  nonlinear  Schrbdinger  equation  becomes  multi-stable.  This 
implies  that  more  than  one  amplitude  proGle  and  speed  of  propagation  of  a  singular  sol¬ 
iton  may  exist  for  the  same  amount  of  the  total  power  carried  by  the  soliton.  The 
existence  of  multi-stable  solitons  is  related  to  the  properties  of  nonlinear  susceptibility 
as  a  function  of  the  intensity  of  light.  For  example,  the  multistable  soliton  waves  can 
not  be  observed  in  a  Kerr-like  nonlinear  medium;  they  may  exist  only  if  the  nonlinear 
component  of  the  susceptibility  as  function  of  intensity  is  either  changing  its  sign  or  its 
derivative  has  a  sufficiently  sharp  peak  (e.g.  is  a  step- like  function). 

The  soliton  bistability  can  result  in  such  effects  as  bistable  (or  multistable,  in  gen¬ 
eral)  self-trapping  of  light  in  media  with  nonlinear  refractive  index  as  well  as  bistable 
propagation  of  short  soliton  pulses  in  nonlinear  optical  fibers  since  both  of  them  may  be 
described  by  the  same  nonlinear  equation.  Both  of  these  effects  may  be  viewed  as  an 
ultimate  manifestation  of  multistable  wave  propagation  since  they  are  based  on  the 
simplest  possible  propagation  configuration.  They  may  also  provide  new  opportunities 
in  the  field  of  optical  bistability!*!.  Indeed,  for  example,  a  bistable  self-trapping  of  light 
provides  a  potential  for  optical  bistable  device  entirely  free  either  from  any  cavity  or 
Fabry-Perot  resonators!*!,  single  nonlinear  interfaces!^!  or  nonlinear  waveguides  formed 
by  the  nonlinear  interfaces!*!,  retroreflection  self-action  effects^,  four-wave  mixing!*!, 
etc.  On  the  other  hand,  since  the  propagation  of  singular  pulses  in  a  homogeneous  non¬ 
linear  medium!*'^!  and  in  nonlinear  fiber  waveguides!*!  is  also  governed  by  a  nonlinear 
Schrbdinger  equation,  these  soliton  pulses  in  the  system  with  an  appropriate  nonlinear¬ 
ity  may  provide  the  first  (to  the  best  of  our  knowledge)  known  opportunity  to  attain  a 
temporal  (or  dynamic)  bistability  as  opposed  to  all  known  kinds  of  optical  bistability 
which  were  so  far  formulated  in  terms  of  steady-state  regimes.  The  very  notion  of 
steady-state  optical  bistability  comes  into  the  inevitable  contradiction  with  the  applica¬ 
tions  most  of  which  assume  fast  pulse  regime  of  operations.  When  exploited  in  a 
dynamic  regime,  such  effects  still  demonstrate  hysteretic  behavior  which,  however,  can 
hardly  be  identified  with  the  original  "adiabatic”,  steady-state  hysteresis.  The  dynamic 
hysteresis  is  more  strongly  affected  by  the  relaxation  processes  rather  than  by  steady- 
state  bistable  states  especially  when  the  total  switching  cycle  has  the  duration  time  of 
the  same  order  as  relaxation  times.  The  truly  dynamic  (or  temporal)  bistability  dis¬ 
cussed  in  this  paper  is  based  on  bistable  pube  shapes  (as  well  as  on  bbtable  duration  of 
the  pubes)  and  offers  a  way  to  resolve  thb  contradiction. 

We  consider  the  generalized  nonlinear  Schrbdinger  equation  for  the  complex  ampli¬ 
tude  of  field  E  in  the  form 

2i5E/dz  +  d*E/ax*  +  Ef(|  E|  *)  =  0  (1) 

where  f(|E|^)  b  an  arbitrary  function  of  the  intensity  |E|^  with  f(0)  =  0.  When 
f(|  E|  ^)  =  o|  E|  ^(o  =  const),  (1)  b  so  called  cubic  nonlinear  Schrbdinger  equation!^’*"**! 

\  i  AUi>4:.  i* 
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(which  corresponds  to  Kerr-nonlioearity  in  optical  propagation).  In  the  case  of  two- 
dimensional  self-trappingl*‘l,  z  is  a  normalized  axis  of  the  soliton  propagation  and  x  is 
a  normalized  transversal  axis  (both  of  them  are  dimensionless  and  correspond  to  the 
real  coordinate  z  and  x  multiplied  by  the  wave  number  k  =  um/c),  whereas  in  the  case 
of  pulse  propagation  (along  the  axis  z),  x  is  interpreted  as  t  -  z/v,  v  =  du/dk  is  a 
group  velocity  in  the  linear  case,  with  the  initial  equation  of  propagation  beingl^^l 

2ikdE/az  +  (dv/du)v-2a2E/dx*  +  k^Efd  E|  *)  =  0  ,  (1‘ ) 

and  can  readily  be  transformed  into  Eq  (1)  by  proper  scaling.  In  both  of  the  cases  f  is 
proportional  to  the  normalized  nonlinear  (i.e.  intensity  dependent)  component  of 
the  dielectric  constant  e  of  medium.  The  nonlinear  Schrodinger  equation  b  obtained 
from  the  Maxwell  equations  in  a  conventional  approximation  of  slowly  varying  envelop 

(i.e.  «  — r)  which  implies  either  small  (quasi-optical)  diffraction  (Eq  (1)J  or 

dz^  dx^ 

relatively  small  dbpersion  [Eq  (!')]. 

The  stationary  solutions  (in  particular,  singular  solitons)  of  Eq  (1)  have  nonvarying 
intensity  profile,  d|E|Vdz~0«  such  solutions  are  written  as 

E(x,z)  =  u(x)exp(j2/2+i^),  where  0  =  const  and  6  b  (unknown)  real  speed  (or  propar 
gational  constant)  of  the  soliton.  Thus,  the  equation  for  the  real  amplitude  tt(x)  b 

d^u/dx*  +  u(f(u2)  -  «)  =  0  (2) 

whose  soliton  solution  must  satbfy  the  condition  u-^^Oas  |x|  -♦ooin  order  for  the 

00 

total  power  P  =  /  u^dx  to  be  limited.  Thb  provides  for  the  first  integral  of  Eq  (2)  in 

-oo 

the  form 

U 

(du/dx)*  =2  /  u|5-  f(u*))du  (3) 

0 

integration  of  which  gives  the  soliton  amplitude  profile  u(x)  for  each  particular  6  and 
f(u^).  In  order  to  evaluate  a  total  power  P,  however,  one  needs  not  to  know  an  explicit 
form  u(x).  Indeed,  by  making  use  of  Eq.  (3)  and  introducing  I  =  u^  =  ]  E|  one  shows 
that 

U<) 

p(«)=  /  di/>/n^,  (4) 

0 

where 

I 

F(I)  =  !-«/  f(l)dl ,  (F(0)=0)  (5) 

0 

i.e.  P  b  determined  immediately  by  f(I)  and  6.  In  Eq  (4),  I,„( £)  b  a  peak  intensity  of 
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the  soliton;  it  is  defined  as  a  minimal  positive  root  of  the  equation  F(I)  =  6.  The 
multistability  of  a  singular  soliton  b  realized  when  the  function  j(P)  implicitly  deter¬ 
mined  by  Eq  (4)  becomes  multivalued. 

It  b  readily  shown  that  a  Kerr-like  nonlinearity  (foci)  results  only  in  a  one-valued 
singular  soliton  (with  see  Fig.  1,  curve  S;  the  same  b  valid  also  for  any  other 

nonlinearity  with  f°cP  where  /i>0  (but  ft  ^2).  The  nonlinearity  f«I^  plays  a  special  role 
in  two-dimensional  propagation  in  the  sense  that  in  thb  case  the  total  energy  carried 
by  any  singular  soliton  b  the  same  regardless  of  its  spatial  profile  and  preparation  con¬ 
stant.  Indeed,  for  f  =  I^Io  where  I,  =  const,  the  intensity  profile  I(x)  and  propagation 
constant  6  are  defined!'^!  from  Eq  (3): 

I(x)  =  I„/co.h(2I„«/lov^) ;  S  =  (6) 

where  maximal  intensity  of  soliton  b  an  arbitrary  constant;  the  total  power  b 
P  =  x/2lo.  One  may  note  from  Eq  (4)  that  in  the  general  case  of  arbitrary  f(l),  a  con¬ 
stant  S  may  be  viewed  as  a  first  integral  (“energy”)  of  some  system  with  a  potential 
F(I),  Eq  (5).  The  motion  of  thb  system  in  some  p  domain  can  then  be  described  by  the 
equation 

d^I/dp*  +  8d{F(I)l/dI  =  0  (7) 

where  if  p  b  interpreted  as  a  “time”,  P(j)  b  a  total  “period”  of  oscillation  of  the  sys¬ 
tem  for  any  given  “energy”  of  excitation  6.  Particularly,  one  may  see  that  the  case 
f«I^  (and  therefore,  FkI^)  corresponds  to  a  “linear  oscillator”,  with  the  “period”  of 
oscillation  P  independent  on  its  “energy”  6,  i.e.  dP/d^  =  0  as  suggested  above. 

In  order  to  demonstrate  the  exbtence  of  the  countable  set  of  states  of  the  singular 
soliton  (with  more  than  one  state)  we  consider  first  the  step- nonlinearity: 


f(I)  =  0,  if  I  <  Iq  ,  and  f  =  A  ,  if  I  >  1q 
where  1q  and  A  are  some  positive  constants.  Substituting  (8)  into  (4)  one  gets 
prn  ~  ^  ^  I  ^  +  «csin\/^  V  f,—  6 


(8) 

(0) 


The  function  ^  vs  P  determined  by  (0)  b  a  two-valued  function  (Fig.  1,  curve  1)  for  any 
P  >  P„i  with  /9(Pcr)  ^0.21.  The  further  example  b  given  by  the  non¬ 

linearity 


f(I)  =  0  ,  if  I  <  lo,  and  ni)  =  Ml  -  loW  .  J  >  *0  •  (W) 

In  contrast  with  (8),  f(I)  b  now  a  continuous  function,  whereas  its  derivative  df/dl  b 
still  dbcontinuous.  The  total  power  (4)  now  b 
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arc  cos 


); 


(11) 


which  essentially  represents  the  same  kind  of  behavior  as  (0),  i.e.  provides  two-valued 
solution  ^(P)  for  any  P  >  P„2  w4.28Iq/>/A  with  ^{P„)  «0.26.  In  these  cases,  the 
nontrivial  branches  of  the  function  P(j)  tend  to  infinity  as  £  0  and  6  -*A  (note  that 

the  third,  “trivial”,  branch  with  £=0,  P-arbitrary,  corresponds  to  a  nontrapped  beam 
with  <  lo)-  This  suggests  a  bistability  without  hysteresis  and  is  due  to  the  fact  that 
nonlinearity  f(I)  differs  from  zero  only  for  some  finite  I  >  Iq.  The  same  kind  of  soliton 
bistability  is  exhibited  by  the  system,  if  either  (i)  f  (0)  <  0  but  f(I)  becomes  positive  at 
some  I,  e.g.  when  f  =  -a,I  +  aji*  -  aji*  where  a|,a2,a3  >  0  and  Qaiaj  <  2a|,  or  (ii) 
f(I)  >  0  in  the  vicinity  of  I  =  0  but  f(I)  =  o{I*);  e.g.  f(I)  =  a, I®  -  a2l^a|,a2  >  0)  or 
f(I)  =  a,I®/(l  +  I®/!®).  (ai,Io  >  0).  The  latter  nonlinearity  may  result  from  the  three- 
photon  resonant  absorption  of  light  by  two-level  systems  with  saturation. 

In  order  to  attain  truly  hysteretic  bistable  behavior  (i.e.  that  characterized  by  the 
S-shape  steady-state  curves  (see  e.g.  curves  2  and  3  at  Fig.  1)  which  causes  both  “on” 
and  “off”  jumps  between  different  branches  of  the  curve],  the  function  f(I)  must  be 
positive  at  least  in  some  range  0  <  I  <  I|  and  have  a  distinct  peak  of  its  first  deriva¬ 
tive  r  (I)  in  this  range.  The  existence  of  hysteretic  jumps  is  secured  if  dS/dP  =  oo  (or 
dP/dS  =  0)  for  two  (or  more)  discrete  values  of  P  (or  6),  where  dP/dS  is  found  from  (4) 
as 


dP  ^  1  .F(d®F/di®i ,  di 

dS  26  y  (dF/dl)*  v'tf-Fd) 


(12) 


A  derivative  "^(^j  is  strongly  affected  by  F' '  (I)  and  therefore  by  f  (I);  bistability  may 
dd 

exist  if  f'(O)  >  0,  and  if  at  some  point  1=1,  there  is  P'(i)  =0  and  r(i)  >  r(0).  As 
an  example  of  such  a  function,  consider 


f  =  aj  +  agl*  -  a3l® 


(13) 


where  ai,a2,a3  >  0.  S-shape  behavior  of  P(^)  (Fig.l,  curve  3)  is  possible  if  the  condition 
is  satisfied 


a,a3/a|  <  S,,  =  0(1)  (14) 

where  S„  is  some  critical  quantity;  the  rough  estimate  gives  S„  ~0.1-0.2.  In  general 
case,  the  critical  situation  (when  the  curve  P(6)  at  some  point  6  =  6„  has 

=  0,  see,  e.g.  Fig.l,  curve  4)  corresponds  to  the  conditions 

dd  dir 

dP/d6„  =0  and  2  (d*F/dl®)  F  =  (dF/dI„)*  (15) 

(where  I^,  is  the  minimal  solution  of  the  equation  6„  =  F(I„)]>  which  determines  both 


6„  and  the  required  parameters  of  the  function  F(I)  [and  therefore,  f(I)].  In  the  case 
when  f(I)  =  0(1*)  at  1=0,  the  function  6(P)  forms  a  hysteresis  if  d*f/dl*  >  0; 
d*f/dl*  >  0,  and  d^f/dl^<0  at  1=0,  e.g.,  f  =  a2l*  +  a3l*-a4l*  (a2,a3,a4  >  O).  In  such  a 
case,  the  lower  (stable)  branch  of  ^P)  corresponds  to  non-trapped  beam  (£=0),  see  Fig. 
1,  curve  2. 

A  stability  of  each  of  possible  solitons  which  correspond  to  the  same  total  power  P 
is  an  important  issue.  The  detailed  small  perturbation  analysis  of  the  spatial  stabil¬ 
ity  of  multistable  solitons  in  the  case  of  step-nonlinearity  (8)  shows  that  the  lower 
branch  of  curve  1,  Fig.  1  corresponds  to  the  unstable  solitons  and  the  upper-  to  the 
stable  ones;  the  trivial  solution  (£  =  0)  is  stable  for  any  P,  This  suggests  a  general  cri¬ 
terion  for  an  arbitrary  f(l),  and  therefore  ^P):  the  stable  solitons  are  those  for  which 
dd/dP  >  0  and  vice  versa  (see  Fig.  1,  curves  1-3).  In  a  further  study,  it  is  of  consider¬ 
able  interest  to  study  a  "collision”  of  two  solitons  that  belong  to  upper  and  lower 
branches  of  the  curve  d(P). 

The  bistable  solitons  may  exist  also  in  the  case  of  three-dimensional  propagation. 
Stationary  self-trapping  of  a  cylindrical  beam,  for  instance,  is  governed  by  the  "non¬ 
linear  Bessel”  equation  [instead  of  (2)|: 

d*u/dr*  +  (1/r)  (du/dr)  +  u[f(u*)  -  d)  =  0  (16) 

where  r  is  a  radial  coordinate  in  the  plane  normal  to  z  axis.  For  cylindrical  beams,  a 
Kerr-nonlinearity,  f«I,  plays  the  same  role  as  f<xl*  in  two  dimensional  case;  for  such  a 
nonlinearity,  the  total  power  of  the  beam  dora  not  depend  of  its  size  or  its  peak  inten¬ 
sity  1*1.  Therefore,  in  order  to  attain  a  nonhysteretic  bistable  soliton  propagation  of  the 
kind  depicted  by  curve  1,  Fig.  1,  the  lower  required  degree  of  nonlinearity  at  I  —»  0  is 
f«I^  [with  f  attaining  some  maximum  or  saturation  when  I  increases,  e.g. 
f  =  aI*/(l+F/I^)|.  Such  a  nonlinearity  can  be  originated  e.g.  by  the  two-photon 
resonant  absorptionl^^l.  Hysteretic  characteristic  curve  d(P)similar  to  curve  2,  Fig.  1, 
can  be  provided  now  by  the  nonlinearity  f(I)  =  a|l  +  a2l*  -  a3l*  (a|,a2,a3  >  0),  with 
the  critical  condition  in  the  same  form  as  Eq  (14)  [but  with  different  S„  =  0(1)]. 

In  conclusion,  we  demonstrated  an  existence  of  multi-stable  solitons  of  generalized 
nonlinear  Schrodinger  equation.  In  order  for  those  solutions  to  exist,  the  nonlinearity 
must  satisfy  some  special  conditions,  e.g.  its  dependence  on  the  light  intensity  must 
have  a  range  where  it  increases  sufficiently  sharply.  In  nonlinear  optics,  these  solitons 
may  manifest  themselves  either  as  singular  pulses  (e.g.  in  nonlinear  fibers)  or  self- 
trapped  channeb  (both  in  two-  and  three-dimensional  cases).  Bbtable  solitons  present 
the  ultimate  case  of  multbtable  wave  propagation  and  may  find  an  application  to  the 
dynamic  (temporal)  optical  bistability  and  bbtable  resonator-free  self-trapping  of  light. 

Thb  work  was  supported  by  the  US  Air  Force  Office  of  Scientific  Research. 
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Captions  to  Figure 


Fig.  1.  A  propagation  constant  S  vs  the  total  power  P  carried  by  the  soliton.  Curves 
1*5  corresponds  to  various  functions  of  nonlinearity:  1  >  step-function  Eq  (8);  2 
-  f  =  a2l^+a3p-a4p(a2,a3,a4>0):  3  -  Eq  (13)  with  aja,  <  a|S„;  4  -  Eq  (13) 
with  ajaj  =  5  -  Kerr-noniinearity,  f  1.  The  broken  lines  at  curves  1-3 

correspond  to  the  unstable  solitons.  In  the  insertion,  the  intensity  profiles  I(x) 
are  depicted  of  solitons  that  carry  the  same  power  but  correspond  to  different 
branches  of  function  ^P)  -  upper  branch  (U)  and  lower  branch  (L). 
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Abstract 

It  is  demonstrated  that  a  nonlinear  Schrodinger  equation  with  certain  nonlineari¬ 
ties  allows  for  an  existence  of  multi-stable  single  solitons  (i.e.,  single  solitons  with  the 
same  carried  power  but  different  propagation  parameters).  In  nonlinear  optics,  these 
solitons  may  exist  either  in  the  form  of  short  bistable  pukes,  or  bistable  self-trapping 
(both  two-  and  three-dimensional). 
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1.  Introduction 

The  most  known  optical  bistability  principles  and  devices  are  based  on  Fabry- 
Perot  resonators  filled  by  a  nonlinear  material.  The  entirely  different  principle  of 
switching  effect  has  first  been  proposed  1^1  and  subsequently  observed  in  the  experi¬ 
ment  1^1  which  does  not  exploit  any  resonant  effect  and  is  based  on  the  reflection  of  light 
at  the  single  nonlinear  interface  formed  by  semi-infinite  linear  and  nonlinear  media.  In 
the  recent  work  W  (with  the  earlier  proposals  1^1)  this  principle  has  been  extended  to  the 
so-called  nonlinear  waveguides  formed  by  two  nonlinear  interfaces  (i.e.  sandwich  struc¬ 
ture  with  a  linear  layer  between  two  semi-infinite  nonlinear  layers  or  vice  versa  1^1). 
Although  the  nonlinear  waveguides  suggest  new  interesting  opportunities,  they  are  still 
substantially  based  on  the  nonlinear  interface  principle. 

In  this  paper,  the  feasibility  of  a  fundamentally  novel  bistable  and  switching  opti¬ 
cal  effect  is  demonstrated  which  consists  in  a  multivalued  self-trapping  of  light  in  an 
inttrface-free  nonlinear  medium.  For  the  certain  class  of  nonlinearities  the  self-trapped 
beam  of  light  occurs  to  have  more  than  one  possible  propagation  constant  (and  respec¬ 
tively  more  than  one  possible  intensity  profile  in  its  cross-section)  for  the  same  total 
power  carried  by  the  beam. 

Self- trapping  and  self- focusing  of  light  in  nonlinear  media  may  be  described  by  a 
parabolic  nonlinear  partial  differential  equation  (the  so-called  nonlinear  Schrodinger 
equation)  whose  single  solutions  (solitons  t**))  correspond  to  self-trapped  channels 
We  show  that  the  existence  of  multi-soliton  solutions  of  the  (generalized)  Schrodinger 
equation  is  related  to  the  properties  of  nonlinear  susceptibility  as  a  function  of  the 
intensity  of  light.  For  example,  the  multistable  soliton  waves  can  not  be  observed  in  a 
Kerr-like  nonlinear  medium;  they  may  exist  only  if  the  nonlinear  component  of  the  sus¬ 
ceptibility  as  function  of  intensity  is  either  changing  its  sign  or  its  derivative  has  a 
sufficiently  sharp  peak  (e.g.  is  a  step-like  function).  This  may  be  e.g.  due  to  such 


mechaDisms  as  either  light-induced  phase  transition  in  a  material  or  multi-photon 
saturation  of  resonant  levels.  This  special  requirement  may  probably  be  a  reason  why 
the  multi-stable  soliton  effects  has  not  been  dbcussed  earlier. 

Furthermore,  since  the  propagation  of  single  pulses  in  a  homogeneous  nonlinear 
mediuml^’**!  and  in  nonlinear  fiber  waveguidesl’^1  is  also  governed  by  a  nonlinear 
Schrodinger  equation,  these  soliton  pulses  in  the  system  with  an  appropriate  nonlinear¬ 
ity  may  provide  the  first  (to  the  best  of  our  knowledge)  known  opportunity  to  attain  a 
temporal  (or  dynamic)  bistability  as  oppmed  to  all  known  kinds  of  optical  bistability 
which  were  so  far  formulated  in  terms  of  steady-state  regimes.  The  very  notion  of 
steady-state  optical  bistability  comes  into  the  inevitable  contradiction  with  the  applica¬ 
tions  most  of  which  assume  fast  pulse  regime  of  operations.  When  exploited  in  a 
dynamic  regime,  such  effects  still  demonstrate  hysteretic  behavior  which,  however,  can 
hardly  be  identified  with  the  original  “adiabatic",  steady-state  bysteresb.  The  dynamic 
hysteresb  is  more  strongly  affected  by  the  relaxation  processes  rather  than  by  steady- 
state  bistable  states  especially  when  the  total  switching  cycle  has  the  duration  time  of 
the  same  order  as  relaxation  times.  The  truly  dynamic  (or  temporal)  bbtability  db¬ 
cussed  in  this  paper  b  based  on  bbtable  pube  shapes  (as  well  as  on  bbtable  duration  of 
the  pulses)  and  offers  a  way  to  resolve  thb  contradiction. 

Both  of  these  effects  (multbtable  self-trapping  of  light  and  multbtable  soliton 
pubes)  may  be  viewed  as  an  ultimate  manifestation  of  multbtable  wave  propagation 
since  they  are  based  on  the  simplest  possible  propagation  configuration.  They  may  also 
provide  new  opportunities  in  the  field  of  optical  bbtability. 


2.  Wave  Equation. 

We  assume  that  the  EM  field  E(f)e"'^^‘)  propagates  in  a  lossless  medium  having 
intensity-dependent  susceptibility  £,  such  that 

+  A€NL(|  E|  2)  (1) 

and  introduce  f(|  E|  2)  =  where  f  can  be  an  arbitrary  function  of  the  intensity 

I  E|  2  with  f(0)  =  0.  The  case  of  f  «  |  Ej  *  corresponds  to  the  Kerr-nonlinearity.  We 
assume  also  that  the  axis  of  the  propagation  is  z  (i.e.  T  =  kCj,  where  k*  =  and 

introduce  dimensionless  coordinates  z  =  kz,  x  =  kx,  y  =  ky.  Then  from  the  Maxwell 
equation,  in  the  conventional  approximation  of  a  slowly  varying  envelope 
(d'^/dz^  <<  AjE,  where  Ej^  is  a  Laplacian  operator  in  a  plane  normal  to  the  z  axis), 
one  readily  gets  the  (generalized)  nonlinear  Schrbdinger  equation  governing  the  non¬ 
linear  wave  propagation; 

2iaE/az  +  A^E  +  Ef(l  E]  2)  =  0  (2) 

In  the  two-dimensional  case  (e.g.  with  djdy  -  0),  this  equation  is  reduced  to  the  form 

2iaE/dz  +  a2E/ax2  -I-  Ef(|  E|  2)  =  0  (3) 

In  the  case  of  the  one-dimensional  pulse  propagation  along  the  Z]  axb  in  the  slightly 
dispersive  medium  with  a  nonlinearity  f|(jE|2),  the  equation  of  the  propagation 

isl«.»».»2l 

2i5E/dz,  +  (dv/dw)v-2a2E/de2  +  kEf,(|  E|  *)  =  0  (4) 

where  ^  =  t-Zj/v;  v  =  da;/dk  is  the  group  velocity  of  linear  propagation.  Eq.  (4)  can 
readily  be  transformed  into  Eq.  (3)  by  proper  scaling,  e.g.  by  assuming 

Z|  =  z/k2(dv/dw);  ^  =  x/kv;  fj  =  f k(dv/dfa;).  (5) 

The  slow  varying  envelope  approximation  implies  either  small  diffraction  (Eqs.  (2),  (3)] 
or  small  dispersion  [Eq.  (4)]. 
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The  stationary  solutions  (in  particular  single  solitons)  of  Eqs.  (2-4)  have  nonvary¬ 
ing  intensity  profile,  d|  E|  ^/dz  =0,  i.e.  such  solutions  are  written  in  the  form 

E(x,y,z)  =  u(x,y)xp(Wz/2  +  i^)  (6) 

where  u  =  |E|  is  a  real  amplitude  of  the  field,  ^  =  count  is  a  real  phase  and  6  is  (unk¬ 
nown)  real  speed  (or  propagation  constant)  of  the  soliton. 

3.  Bistable  Soliton  Pulses  and  Two-Dimensional  Self-Trapping 

Both  two-dimensional  self-trapping  and  one-dimensional  soliton  pulses  (e.g.  in 
fibers  1*^1)  are  governed  by  the  same  Eq.  (3).  Substituting  Eq.  (6)  into  Eq.  (3)  and  bear¬ 
ing  in  mind  that  in  this  case  u  =  u(x),  one  gets  the  equation  for  the  real  amplitude  u: 

d^u/dx*  +  *»(f(u^)  -  ^  =  0  (7) 

whose  soliton  solution  must  satisfy  the  condition  u-»0as  |x|  -+ooin  order  for  the 

00 

total  power  P  =  /  u*dx  to  be  limited.  This  provides  for  the  first  integral  of  Eq  (7)  in 

-00 

the  form 

U 

(du/dx)*  =  2  /  u(«  -  f(u2)ldu  (8) 

0 

integration  of  which  gives  the  soliton  amplitude  profile  u(x)  for  each  particular  6  and 
f(u^).  In  order  to  evaluate  a  total  power  P,  however,  one  needs  not  to  know  an  explicit 
form  u(x).  Indeed,  by  making  use  of  Eq.  (8)  and  introducing  I  =  u^  =  |  E|  one  shows 
that 

U<) 

p(«)  =  /  <o/v/ni5) ,  («) 

0 


where 
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I 

F(I)  =  r»/  f(I)dl  ,  (F(0)=0)  (10) 

0 

i.e.  P  is  determined  immediately  by  f(I)  and  S.  In  Eq  (0),  I„(j)  is  a  peak  intensity  of 
the  soliton;  it  is  defined  as  a  minimal  positive  root  of  the  equation  F(I)  =  6.  The 
multistability  of  a  single  soliton  is  realized  when  the  function  j(P)  implicitly  determined 
by  Eq  (0)  becomes  multivalued. 

It  is  readily  shown  that  a  Kerr-like  nonlinearity  (f°cl),  Fig.  2,  curve  1,  results  only 
in  a  one-valued  single  soliton  (with  see  Fig.  1,  curve  5;  the  same  is  valid  also  for 

any  other  nonlinearity  with  f«I^  where  /<>0  (but  /<^2).  The  nonlinearity  f“cl*  Fig.  2, 
curve  2,  plays  a  special  role  in  the  two-dimensional  propagation  in  the  sense  that  in  this 
case  the  total  energy  carried  by  any  single  soliton  is  the  same  regardless  of  its  spatial 
profile  and  propagation  constant.  Indeed,  for  f  =  I*/Ip  where  I,,  =  const,  the  intensity 
profile  I(x)  and  propagation  constant  S  are  defined!'^)  from  Eq  (8); 

I(x)  =  I„/c<»h(2l„x/Iov/3) ;  6  =  I’/3Io*  (11) 

where  maximal  intensity  of  soliton  is  an  arbitrary  constant;  the  total  power  is 
P  =  7r/2lo.  One  may  note  from  Eq  (fi)  that  in  the  general  case  of  arbitrary  f(l),  a  con¬ 
stant  6  may  be  viewed  as  a  first  integral  ("energy”)  of  some  system  with  a  potential 
F(I),  Eq  (10).  The  motion  of  this  system  in  some  p  domain  can  then  be  described  by 
the  equation 

d*I/<ip*  +  8<l(F(I)l/dI  =  0  (12) 

where  if  p  is  interpreted  as  a  "time”,  P(j)  is  a  total  "period”  of  oscillation  of  the  sys¬ 
tem  for  any  given  "energy”  of  excitation  b.  Particularly,  one  may  see  that  the  case 
fecp  (and  therefore,  F°cl^)  corresponds  to  a  "linear  oscillator”,  with  the  "period”  of 
oscillation  P  independent  on  its  "energy”  b,  i.e.  dP/db  =  0  as  suggested  above. 


f(I)  =  a,I*  -  a2l^(a,,a2  >  0)  or  f(I)  =  a,lV(l  +  I*/**).  >  0).  F»8-  2,  curve  6.  The 

latter  nonlinearity  may  result  from  the  three-photon  resonant  absorption  of  light  by 
two-level  systems  with  saturation. 


In  order  to  attain  truly  hysteretic  bistable  behavior  [i.e.  that  characterized  by  the 
S-shape  steady-state  curves  (see  e.g.  curves  2  and  3  at  Fig.  1)  which  causes  both  “on” 
and  “off”  jumps  between  different  branches  of  the  curve],  the  function  f(I)  must  be 
positive  at  least  in  some  range  0  <  I  <  I|  and  have  a  distinct  peak  of  its  first  deriva¬ 
tive  df/dl  in  this  range.  The  existence  of  hysteretic  jumps  is  secured  if  d^/dP  =  oo  (or 
dP/d5  =  0)  for  two  (or  more)  discrete  values  of  P  (or  5),  where  dP/d6  is  found  from  (9) 
as 


iP  -  1  n,  ;F(,l»F/dP)  ■  dl 
2S  y  {dF/dlf  VPm 


(17) 


A  derivative  -77(5)  is  strongly  affected  by  d*F/dI®  and  therefore  by  df/dl;  bistability 
do 

may  exist  if  df/dl  >  0,  and  if  at  some  point  1  =  1,  there  is  d*f/dl*  =  0  and 
df/dl  >  df(0)/dl.  As  an  example  of  such  a  function,  consider  (Fig.  2,  curve  7) 


f  =  ail  +  aji*  -  ajl^ 


(18) 


where  ai,a2,a3  >  0.  S-shape  behavior  of  ^P)  (Fig.l,  curve  3)  is  possible  if  the  condition 
is  satisfied 


aiaj/al  <  S„  =  0(1)  (19) 

where  is  some  critical  quantity;  the  rough  estimate  gives  S„  ~  0.1-0.2.  In  general 
case,  the  critical  situation  [when  the  curve  P(j)  at  some  point  6  =  has 

=  0,  see,  e.g.  Fig.l,  curve  4]  corresponds  to  the  conditions 

dd  dr 

dP/d6„  =  0  and  2  (d*F/dI*)  F  =  (dF/dI„)*  (20) 

[where  I„  is  the  minimal  solution  of  the  equation  6„  =  F(I„)],  which  determines  both 


6^^  and  the  required  parameters  of  the  function  F(l)  [and  therefore,  f(I)|.  In  the  case 
when  f(I)  =  0(1*)  at  1=0,  the  function  5(P)  forms  a  hysteresis  if  d*f/dl*  >  0, 
d®f/dl*  >  0,  and  dV/dl*  <  0  at  1=0,  e.g.,  f  =  aji*  +  aji*  -  a4l^,  (a2,a3,a4  >  0),  sec 
Fig.  1,  curve  2.  In  such  a  case,  the  lower  (stable)  branch  of  ^P)  corresponds  to  non- 
trapped  beam  (^=0). 

To  analyze  the  conditions  (20)  on  the  general  case  of  an  arbitrary  function  f(I)  is 
not  an  easy  task  since  the  dependence  P(6)  is  implicitly  determined  by  the  integral  (0) 
which  is  not  evaluated  in  the  analytic  form  for  an  arbitrary  f(I)  [and  therefore  F(I),  Eq. 
(10)].  Therefore,  for  the  practical  purpose,  it  is  important  to  have  a  good  analytic 
approximation  of  the  function  P(j)  at  least  in  the  range  0  <  I  <  Ig  where  Ig  is  the 
point  at  which  d*F(Ig)/dI*  =  0  [one  may  see  from  Eq.  (20)  that  the  critical  point  I  =  I^, 
is  located  in  this  range,  i.e.  I^,  <  Igj.  For  such  a  purpose,  it  is  more  convenient  to 
operate  with  I,„  rather  than  6  [remember  that  6  =  F(I„)j.  Then,  if  F(In,)  is  positive  and 
monotonically  increasing  in  (0,lg),  a  good  approximation  is  given  by  the  formula 


where  F„  =  F(I,„)  =  ^;  F^  =  dF(I„,)/dI;  Fq  =  dF(0)/dI,  and  C  >  0  is  some  constant 
of  the  order  of  unity  which  is  determined  by  the  type  of  the  function  F(I).  For  the  cal¬ 
culations  related  to  the  conditions  (20),  this  constant  is  not  important,  since  the  first  of 
Eqs.  (20)  is  to  be  replaced  now  by  the  condition  dP/dle^  =  0  which  makes  insignificant 
any  scaling  of  P(In,). 

4.  Stability  of  Two-Dimensional  Btetable  Solitons. 

A  stability  of  each  of  possible  solitons  which  correspond  to  the  same  total  power  P 
is  an  important  issue.  Assuming  that  the  amplitude  profile  u,(x)  of  the  particular  soil- 
ton  [determined  by  Eq  (8)]  is  known,  we  represent  a  solution  of  Eq  (2)  in  the  form  of 
perturbed  soliton  solution  (6): 


(22) 


E(x,z)  =  (u,(x)  +  Au(x,z))exp{Wz/2  +  i^(x,z)l 

where  Au,  d^jdz  and  d^/dx  are  small  real  perturbations  and  assume  the  factorization 
of  the  perturbation  in  the  form 

Au  =  w(z)(u'(x)p(x)l;  ^  =  V<z)f(x),  (23) 

where  w,  p,  Vi  ^’^d  (  are  some  unknown  real  functions  of  a  single  variable.  Substitut¬ 
ing  Eqs.  (22)  and  (23)  into  Eq.  (3),  linearizing  Eq.  (3)  with  respect  to  the  small  pertur¬ 
bations,  equalizing  the  real  and  imaginary  parts  of  the  obtained  equation  separately  to 
zero,  separating  the  terms  that  depend  only  on  x  or  z,  and  making  use  of  Eq.  (7),  one 
finally  gets  the  equations  for  the  unknown  functions  w,  p,  and 

dw/dz  =  1/^*;  d^/dz  =  Xw  (24) 

(p'»iV  =MO’(i  (25) 

(26) 

where  i/  and  X  are  some  unknown  real  constants,  and  prime  denotes  now  a  derivative 
with  respect  to  x.  Eqs.  (25)  and  (26)  can  be  also  reduced  to  a  single  equation  of  the 
fourth  order  for  one  variable,  e.g.  p: 

u?((p'»;*)7(0’rr  =t(».*)’p  (2?) 

where  qr  =  i/X.  The  eigenmodes  of  Eqs.  (25-27)  must  satisfy  a  condition 

pu'  -»0,  -»0,  as  |x|  -►  oo.  (28) 

If  these  eigenmodes  as  well  as  their  respective  eigenvalues  v  and  X  (and  therefore 
7  =  i/X)  are  found,  the  sufficient  requirement  for  the  stability  of  the  respective  soliton 
E,  =  u,(x)exp(ijz/2)  is  that  all  7  must  be  negative.  The  soliton  is  unstable  if  any  of  7 
is  positive  (the  sufficient  condition);  the  case  with  7  =  0  has  to  he  further  investigated 


although  normally  it  corresponds  to  a  stable  soliton.  The  lower  order  eigenmodes  of 
Eqs.  (25-27)  with  ^  =  0  are  readily  found.  These  are; 


(i)  X  =  1/  =  0,  p  =  const  0  (i.e.  Au  =  Cu^(x),  C  <<  1),  ^  =  0.  This  mode 
corresponds  to  the  shift  of  the  spatial  position  of  the  soliton  peak  at  the  x  axis  (with  all 
other  characteristics  of  the  soliton  being  intact). 

(ii)  X  =  1/  =  0,  (  =  const  ^  0,  p  ~  0,  which  corresponds  to  the  shift  of  the  soliton 
phase  0. 

(iii)  X  =  0;  ^  0;  p  =  const  ^  0;  (  =  (i/p)x.  This  corresponds  to  the  shift  of  the 

“angle”  of  the  soliton  propagation. 

(iv)  X  9^  0;  ^'  =  0;  ^  =  const  0;  and 


u. 


dx 


0  U, 


^F(u2) 


(29) 


Thb  eigenmode  corresponds  to  the  small  change  of  the  propagation  constant  of  the  soli¬ 
ton  S  (and  therefore,  the  total  power  carried  by  the  soliton),  with  conservation  of  the 
single  soliton  character  of  the  entire  propagation.  Analytical  solution  for  the  higher 
modes  in  the  case  of  arbitrary  f(I)  is  still  to  be  found.  The  detailed  analysis  in  the  case 
of  step-nonlinearity  (13)  shows  that  the  lower  branch  of  curve  1,  Fig.  1  corresponds  to 
the  unstable  solitons  and  the  upper-one  to  the  stable  ones;  (the  trivial  solution  (^  =  0) 
is  stable  for  any  P].  This  suggests  a  general  criterion  for  an  arbitrary  f(I),  and  therefore 
^P):  the  stable  solitons  are  those  for  which  d^dP  >  0  and  vice  versa  (see  Fig.  1, 
curves  1-3).  Although  this  statement  seems  to  be  intuitively  almost  obvious,  it  is  still 
to  be  proven.  It  is  also  of  considerable  interest  to  study  a  “collision”  of  two  solitons 
that  belong  to  upper  and  lower  branches  of  the  curve  ^P). 


6.  Bistable  Three-Dimensional  Self-Trapping 


The  bistable  solitons  may  exist  also  in  the  case  of  three-dimensional  propagation. 
Stationary  self-trapping  of  a  cylindrical  beam,  for  instance,  is  governed  by  the  “non¬ 
linear  Bessel”  equation  which  follows  from  Eqs.  (2)  and  (6)  when  u(x,y)  has  a  polar 
symmetry: 


d^u/dr*  +  (1/r)  (du/dr)  +  u(f(u*)  -  rf]  =  0  (30) 

where  r  =  \/x*+y*  is  a  radial  coordinate  in  the  plane  normal  to  z  axis.  With  the 
step-nonlinearity  (13)  the  solution  of  Eq  (30)  is 


u(r)  = 


|a  Jo(kr%/A-tf)  ,  when  u  >  Uq  =  . 

|b  Ko(krv^)  ,  when  u  <  Uo  ,  r  >  ro 


r  <  To 


(31) 


where  and  K„  are  lA-order  ordinary  and  modified  Bessel  functions  respectively;  A,  B, 
and  ro  are  some  constants.  By  requiring  continuity  of  u(r)  and  du/dr  at  u  =  Uq,  one 
gets  an  equation  determining  the  radius  of  the  beam  rp  at  the  level  u  =  Uq: 

/a  ^  Ji(kro>/A=r)  _  Ki(krov^) 

V  6  Jo(kroN/^)  Ko(kro>/tf)  ’  '  ^ 

With  ro  known,  A  and  B  are  determined  as  A  =  Uo/Jo(kro\/A-tf);  B  =  Uo/Ko(kro^/j). 
For  cylindrical  beams,  a  Kerr-nonlinearity,  f°cl,  plays  the  same  role  as  f°cl^  in  the  two 
dimensional  case:  for  such  a  nonlinearity,  the  total  power  of  the  beam  does  not  depend 
of  its  size  or  its  peak  intensity^^l.  Therefore,  in  order  to  attain  a  nonhysteretic  bistable 
soliton  propagation  of  the  kind  depicted  by  curve  1,  Fig.  1,  the  lower  required  degree  of 
nonlinearity  at  I  0  is  f«I^  [with  f  attaining  some  maximum  or  saturation  when  I 
increases,  eg.  f  =  QfI*/{l+I*/Io)l  which  resembles  curve  6,  Fig.  2.  Such  a  nonlinearity 
can  be  originated  e.g.  by  the  two-photon  resonant  absorption^^^l  S-shape  hysteretic 
characteristic  curve  ^P)  can  be  provided  now  by  the  nonlinearity 
f(I)  =  ajl  +  aji*  -  a3l®  (a|,ao,a3  >  0)  which  resembles  curve  7,  Fig.  2,  with  the  critical 


3- 


conditioD  in  the  same  form  as  Eq  (19)  [but  with  different  =  0(1)]. 

0.  Conclusion. 

In  conclusion,  we  demonstrated  an  existence  of  multi-stable  solitons  of  generalized 
nonlinear  Schrbdinger  equation.  In  order  for  those  solutions  to  exist,  the  nonlinearity 
must  satisfy  some  special  conditions,  e.g.  its  dependence  on  the  light  intensity  roust 
have  a  range  where  it  increases  sufficiently  sharply.  In  nonlinear  optics,  these  solitons 
may  manifest  themselves  either  as  single  pulses  (e.g.  in  nonlinear  fibers)  or  self-trapped 
channels  (both  in  two-  and  three-dimensional  cases).  Bistable  solitons  present  the  ulti¬ 
mate  case  of  multistable  wave  propagation  and  may  find  an  application  to  the  dynamic 
(temporal)  optical  bistability  and  bistable  resonator-free  self-trapping  of  light. 

This  work  was  supported  by  the  US  Air  Force  Office  of  Scientific  Research, 


REFERENCES 


1.  For  the  most  recent  review  on  optical  bistability,  see  E.  Abraham  and  D.  Smith, 
Rep.  Prog.  Phys.  815  (1082).  For  introduction  to  the  field  see  P.  W.  Smith 
and  W.  J.  Tomlinson,  IEEE  Spectrum,  18,  26  (June  1981). 

2.  A.  E.  Kaplan,  JETP  Lett.  114  (1076);  A.  E.  Kaplan,  Sov.  Phys.-JETP  n, 
896  (1977); 

3.  P.  W.  Smith,  J.-P.  Hermann,  W.  J.  Tomlinson,  and  P.  J.  Maloney,  Appl.  Phys. 
Lett.  S5,  846  (1979);  P.  W.  Smith  and  W.  J.  Tomlinson,  IEEE  J.  Quant.  Electr. 
QE-20,  30  (1984). 

4.  H.  Vach,  G.  I.  Stegeman,  C.  T.  Seaton,  and  I.  C.  Khoo,  Opt.  Lett.,  9,  238  (1084). 

5.  P.  W.  Smith  and  W.  J.  Tomlinson,  in  Optical  Bistability,  edited  by  C.  M. 
Bowden,  M.  Ciftan,  and  H.  R.  Robl  (Plenum,  NY,  1981),  463-480;  A.  E.  Kaplan, 
ibid,  p.  447-462. 

6.  R.  Y.  Chiao,  E.  Garmire,  and  C.  H.  Townes,  Phys.  Rev.  Letts.  IS,  479  (1964). 

7.  P.  L.  Kelley,  Phys.  Rev.  Letts.  15,  1005  (1965). 

8.  S.  A.  Akhmanov,  A.  P.  Sukhorukov,  and  R.  V.  Khokhlov,  Usp.  Fir.  Nauk  9S,  19 
(1967)  (Sov.  Phys.  -  Usp.  9S,  609  (1968)). 


9.  See,  e  g.  G.  L.  Lamb,  “Elements  of  Soliton  Theory,”  John  Wiley  Sons,  NY,  1980. 

10.  V.  E.  Zakharov  and  A.  B.  Shabat,  Zh.  Eksp.  Teor.  Fiz.  61,  118  (1971)  [Sov. 
Phys.  -  JETP  S4,  62  (1972)). 

11.  L.  A.  Ostrovskii,  Zh.  Eksp.  Teor.  Fiz.  51,  1189  (1966)  (Sov.  Phys.  -  JETP,  24, 
797  (1967)]. 

12.  A.  Hasegawa  and  F.  Tappert,  Appl.  Phys.  Letts.  23,  143  (1973). 

13.  Note  that  the  solution  (11)  is  not  the  same  as  for  well  known  single  solitons  of  the 
cubic  nonlinear  Schrodinger  equation  with  f(I)«I;  in  the  latter  case 
I(x)«l/cosh^(Ax)  where  A  is  some  constant. 

14.  A.  E.  Kaplan,  "On  the  theory  of  a  self-trapped  waveguide  in  a  medium  with  the 
step-nonlinearity,”  Technical  Digest  of  the  VI  National  Conf.  Nonlinear  Optics 
(USSR),  Minsk,  June  27-July  1,  1972,  p.  44-45  (published  by  the  Moscow  State 
Univ.,  Moscow,  1972,  in  Russian). 

15.  V.  S.  Butylkin,  A.  E.  Kaplan,  and  Yu.  G.  Khronopulo,  Zh.  Eksp.  Teor.  Fiz.  59, 
921  (1970)  (Sov.  Phys.  -  JETP  32,  501  (1971)). 


Captions  to  Figures 


Fig.  1.  A  propagation  constant  S  vs  the  total  power  P  carried  by  the  soliton.  Curves 
1-5  corresponds  to  various  functions  of  nonlinearity:  1  -  step-function  Eq  (13); 
2  -  f  =  ajl*  +  ajl*  -  aj^,  (a2,aj,a4  >  0);  3  -  Eq  (18)  with  aiaj  <  a|S„;  4  -  Eq 
(18)  with  aia^  =  ajSj,;  5  -  Kerr-nonlinearity,  f  «  1.  The  broken  lines  at  curves 
1-3  correspond  to  the  unstable  solitons.  In  the  insertion,  the  intensity  profiles 
I(x)  are  depicted  of  solitons  that  carry  the  same  power  but  correspond  to 
different  branches  of  function  S(P)  -  upper  branch  (U)  and  lower  branch  (L). 

Fig.  2.  Various  functions  of  nonlinearity  f  vs.  the  field  intensity  I  (see  detaik  in  the 
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ABSTRACT 


It  was  recently  predicted  by  us  that  cyclotron  resonance  of  free  electrons  in 
vacuum  and  conduction  electrons  in  semiconductors  may  exhibit  bistable  and  hysteretic 
behavior  which  is  due  to  relativistic  mass-effect  (or  pseudo  relativistic  -  in  semiconduc¬ 
tors).  Based  on  this  prediction,  the  hysteretic  cyclotron  resonance  of  a  trapped  single 
electron  in  vacuum  has  recently  been  experimentally  observed  by  G.  Gabrielse  et  al.  In 
this  paper  we  consider  this  phenomenon  as  an  ultimate  bistability  since  it  is  based  on 
the  most  fundamental  mechanism  of  nonlinearity  (relativistic  mass-effect),  involves  the 
interaction  of  EM  wave  with  the  simplest  single  elementary  particle,  and  demonstrate 
the  first  known  intrinsic  bistability  with  no  macroscopic  optical  feedback.  We  also 
show  that  a  hysteretic  resonance  of  a  single  electron  based  on  relativistic  effects  is  feasi¬ 
ble  also  in  one-dimensional  parabolic  potential  (with  no  magnetic  field  required  to 
attain  a  resonance). 


r- 


Introduction 


Optical  bistability*  is  a  rapidly  expanding  and  promising  field  in  nonlinear  optics 
which  offers  both  new  insights  in  nonlinear  interactions  of  light  with  matter  and  poten¬ 
tials  for  superfast  switching  devices  for  optical  computers  and  optical  signal  processing. 
Therefore,  the  fundamental  physical  problems  related  to  that  phenomenon  have  become 
important  as  well.  One  of  the  most  intriguing  questions  is:  what  is  the  ultimate  physi¬ 
cal  level  of  bistable  interaction  of  light  with  matter?  It  is  feasible  to  realize  (and  possi¬ 
bly,  to  exploit)  the  bistable  interaction  of  the  microscopic  level? 

In  our  recent  work^  it  was  predicted  that  even  a  slight  relativistic  mass-effect  of  a 
single  free  electron  may  result  in  large  nonlinear  effects  such  as  hysteresis  and  bistabil¬ 
ity  in  cyclotron  resonance  under  action  of  EM  wave.  The  predicted  effect  may  be 
regarded  as  the  ultimate  and  fundamental  one  in  many  respects: 

(i)  it  suggests  the  bistable  interaction  of  EM  wave  with  the  single  simplest  micros¬ 
copic  physical  object  -  an  electron 

(ii)  the  nonlinearity  that  makes  the  bistable  interaction  possible  is  based  on  one  of 
the  most  fundamental  physical  effects  -  relativistic  change  of  electron  mass 

(iii)  it  was  the  first  proposed  effect  that  offered  bistability  based  on  the  intrinsic  pro¬ 
perty  of  an  microscopic  object  rather  than  on  macroscopic  optical  feedback  in  a 
nonlinear  medium. 

Very  recently,  based  on  the  prediction^,  the  hysteretic  (bistable)  cyclotron  reso¬ 
nance  of  a  free  electron  has  been  discovered  by  G.  Gabrielse  et  al.®  in  experiment  in 
which  a  simple  electron  has  been  trapped  in  a  Penning  trapp  for  the  period  of  time  as 
long  as  10  months^. 

In  Ref  [2]  it  was  also  suggested  that  the  analogous  effect  (i.e.  hysteretic  cyclotron 
resonance)  can  be  expected  also  in  semiconductors  with  narrow  energy  gap  between 
conduction  and  valence  bands.  This  problem  was  later  addressed  in  Ref  (5).  In  sem¬ 
iconductors,  the  mass-effect  required  for  the  hysteretic  resonance,  is  provided  by  the 


nonparabolicity  of  the  semiconductor  conduction  band  which  causes  a  pseudo- 
relativistic  dependence**’^  of  the  effective  mass  of  conduction  electrons  on  their  momen¬ 
tum  or  energy. 

In  the  both  cases  (i.e.  free-space  electrons  and  conduction  electrons),  the  hysteretic 
resonance  is  attributed  to  the  dependence  of  the  cyclotron  frequency  of  forced  oscilla¬ 
tions  on  the  relativistic  mass  of  the  electron,  and  hence  on  its  momentum  (or  kinetic 
energy).  The  critical  condition  for  bistability  to  occur  is  that  the  nonlinear  (relativistic) 
shift  of  the  cyclotron  frequency  must  be  sufRciently  greater  than  the  linewidth  of  the 
resonance.  To  some  extent,  this  effect  resembles  hysteresis  in  a  classic  nonlinear  oscilla¬ 
tor*  or  in  nonlinear  parametric  systems®. 

From  the  electrodynamic  viewpoint,  one  of  the  most  important  features  of  the  new 
effect  is  that  it  is  based  on  the  intrinsic  properties  of  the  microscopic  components,  not 
on  the  macroscopic  optical  feedback.  This  differs  fundamentally  from  all  conventional 
mechanisms  of  optical  bistability*  which  so  far  have  always  been  based  on  macroscopic 
nonlinear  properties  of  the  media.  Nonlinear  change  in  macroscopic  susceptibility 
under  action  of  the  strong  EM  wave  provides  dramatic  change  in  the  optical  condition 
of  propagation  of  this  wave  under  various  special  circumstances  which,  in  turn,  leads 
again  to  the  change  in  the  susceptibility.  This  so-called  optical  feedback  in  nonlinear 
macro-systems  results  in  the  existence  of  multistable  (in  particular,  bistable)  steady 
states.  No  such  optical  feedback  exists  in  the  hysteretic  electron  resonances.  One  of 
the  consequences  of  this  fact  is  that  those  effects  exibit  also  so  called  cavityless  (or 
resonator-free)  bistability*®.  Recently,  some  new  mechanisms  of  optical  bistability 
based  on  the  intrinsic  bistability  analogous  to  bistable  cyclotron  resonance  (in  the  sense 
that  they  are  attributed  to  the  nonlinear  shift  of  resonant  frequency  of  material)  have 
been  proposed'*’*^  and  experimentally  observed**. 

In  this  paper  we  review  the  theory  of  hysteretic  (bistable)  cyclotron  resonance  of  a 
single  free  electron  (Section  1)  and  conduction  electrons  in  semiconductors  (Section  2). 


We  show  also  that  in  fact  the  presence  of  a  magnetic  field  (which  is  required  to  attain  a 
cyclotron  resonance)  is  not  a  necessary  condition  for  a  hysteretic  effect;  one  may  attain 
a  bistability  based  on  relativistic  mass-effect  even  in  any  one-dimensional  oscillator  hav¬ 
ing  parabolic  potential  well  (Section  3). 

1.  Bistable  (hysteretic)  cyclotron  resonance  of  a  firee  electron 

We  shall  demonstrate^  this  effect  for  the  simplest  case  of  a  single  electron 
immersed  in  a  strong  constant  magnetic  field  Ho  and  interacting  with  an  electromag¬ 
netic  field  Ei„  of  amplitude  E  (E«Ho).  The  plane  EM  wave  propagates  along  the  axis 
z  parallel  to  Ho-  Field  Ho  provides  a  cyclotron  resonance  with  the  unperturbed  fre¬ 
quency 

Wo  =  eHo/moC  (1.1) 

where  e  is  an  electrical  charge  of  the  electron,  mo  is  its  rest  mass,  and  c  is  the  speed  of 
light.  We  also  assume  that  a  small  constant  electric  field  ■g(z)  is  applied  along  axis  z; 
this  field  is  provided  by  some  potential  to  arrange  a  trapping  of  the  electron’®  and  to 
compensate  a  radiation  force  caused  by  the  EM  wave. 

We  treat  this  problem  classically.  The  equation  of  motion  for  the  electron  moving 
with  arbitrary  velocity  v  is 

d(mV)/dt  =  (e/c)^  X  fls  eRi:  +  .  (1-2) 

m  =  mo(l  -  1^1  (1.3) 

where  llj;  is  the  total  magnetic  field  (including  the  EM  wave  component  Hemi  * 
Bs  =  n,  +  Bem),  =  Ejn  +  g(z)e„  is  the  total  electric  field,  and  the  term  Tt 
represents  energy  losses  of  the  electron.  In  the  ultimate  case  in  which  the  losses  are 

caused  by  EM  radiation  of  the  rotating  electron  (and  |  v|  «c)  this  term  can  be  writ- 


Ff  =  (2eV3c2)dV/dt2  .  (1.4) 

In  the  general  case  the  losses  may  be  much  larger  and  caused  by  various  factors,  basi¬ 
cally  by  collisions,  so  that  the  radiation  losses  could  be  neglected.  The  force  is  then 
proportional  to  the  velocity  of  the  electron,  e.g.,  F#  =  -TmoWjV,  where  7  is  the  dimen¬ 
sionless  width  of  cyclotron  resonance.  The  radiation  losses  can  also  be  represented  by 
this  formula,  since  one  can  assume*'*  that  dV/dt^  w  which  yields 

2e^a’o  2 

^  =  7  Teko  «  1  .  (1-5) 

3moC'’  o 

where  r,  =  e^/mQC^  =  2.8  x  10*^  cm  is  an  electron  radius  and  kj  =  oiq/c  is  a  resonance 
wave  number.  Since  for  the  plane  wave 

^EM  =  k  X  Ein/w  (1.6) 

(where  w  is  a  frequency  of  EM  field,  and  k  =  w/c),  Eq  (1.2)  governing  the  interaction  of 
the  “magnetized”  electron  with  an  arbitrary  propagating  plane  wave  Ejn  can  be  written 
in  the  form; 


d(mV)/dt  +  THioWoV  =  e 


E„  +  ^rxB, 


+  e 


77  X  (kxEiJ-f 

U; 


(1.7) 


where  the  term  eVx[lcxEin|/w  is  a  radiation  force  applied  to  the  electron.  Usually  this 
term  is  neglected  (see,  e.g.  Ref.  15),  except  in  Ref.  16  in  which,  however,  losses  and 
other  possible  forces  like  counter-potential  X  were  not  considered.  However,  all  of  these 
interactions  become  important^  when  considering  excited  steady  states  (and  multi¬ 
stability)  of  the  electron  under  the  action  of  the  sufficiently  intense  EM  wave. 

We  introduce  dimensionless  notations; 


V 

C 


P  = 


moc^ 


efilil  .  J_  , 

moc^  •'o 


(1.8) 


All  these  variables  (as  well  as  7)  are  supposed  to  be  very  small  compared  to  unity; 


1^1  ,  i/^l  ,  |p|  ,  Tf  «  1  (1.9) 

Using  these  conditions,  and  making  use  of  definitions  (1.8)  and  of  the  fact  that  under 
the  assumed  geometry  the  vectors  flg*  parallel  to  each  other,  one  can 


rewrite  Eq  (1.7)  in  the  form: 


(1.10) 


where  the  term  |  0\  ^/2  represents  a  small  nonlinearity  which  is  due  to  weak  relativistic 
ma.ss  effect.  We  also  a.ssume  that  the  EM  wave  E!;„  is  circularly  polarized  (which  max¬ 
imizes  the  expected  effect)  and  rotates  in  the  same  direction  as  the  electron,  i.e. 

H  =  /i[e,,sin(ut  -  kz)  +  eyCOs(ujt  -  kz))  .  (1.11) 

The  required  solution  to  Eq  (1.10)  can  then  be  written  in  the  form: 

^(t,z)  =  i5(exSin(wt  -  kz  +  ^>)  +  eyCOs(wt  -  kz  +  <;i))  +  /?,e, .  (112) 

By  virtue  of  conditions  (1.9),  the  unknown  variables  0,  and  <(>  vary  little  in  the  time 
u;"',  which  allows  us  to  use  slow-varying  envelope  approximation,  i.e.  to  neglect  their 
second-order  time  derivatives  and  higher  harmonics,  and  to  write  down  the  set  of  trun¬ 


cated  first-order  equations: 

+  H  cos  <(f 


(1.13) 


+  /^/2  +  psin 


(1.14) 


i^z/^0  =  “  p(z)  +  1«/?COS0  ; 


(1.15) 


Z  =C0^ 


(1.16) 


Here  A  =  (w  -  u)o)/u»0  «  1;  A  is  a  dimensionless  resonant  detuning.  The  steady-state 
solution  (d/dt  =  0)  is  thus  determined  by  the  relationships 


^2  =  ^2|^  +  +  ^2/2)2j 


(1.17) 


(,3,),  =  0  ;  /)(Z,)  = 

(1.18) 

tan  =  -(A  +  0^/2)/'i 

(1.19) 

where  the  subscript  “s”  labels  characteristics  of  the  steady-state  regime.  It  is  seen  from 
Eqs  (1.17)  and  (1.19)  that  the  steady-state  kinetic  energy  does  not  depend  on 
parameters  of  the  force  /j(z).  This  force  canceb  the  radiation  force  of  the  incident  wave 
and  is  small  compared  to  all  the  other  forces  in  the  problem;  e.g.,  in  the  region  of 
interest  p  —  ~  0*.  It  is  obvious,  though,  that  the  characteristics  of  the  transient 

regime,  described  by  Eqs  (1.13-1.16)  depend  on  the  spatial  behavior  of  the  force  p(z). 
For  example,  the  frequency  of  transient  oscillations  along  the  z  axis  should  depend 
upon  gradient  dp/dz  in  the  vicinity  of  z  =  z„  Eq.  (1.18).  However,  the  magnitude  of 
these  oscillations  at  the  onset  of  the  hysteretic  jumps  can  be  assumed  to  be  arbitrarily 
small,  provided  that  the  amplitude  (or  the  frequency)  of  the  incident  wave  is  swept 
sufficiently  slow  to  prevent  the  transient  regime  from  masking  the  hysteretic  behavior 
of  the  steady-state  regime.  Under  the  threshold  conditions 

>  pI  =  (16/3>/3)y  ,  A  <  =  -'Tv/S  ,  (1.20) 

Eq.  (1.17)  yields  a  three-valued  solution  for  The  plot  of  steady-state  kinetic  energy 
0^/2  as  a  function  of  resonant  detuning  A  and  incident  intensity  is  shown  in  Fig.  1. 
At  the  threshold  point  this  value  is  -  27/\/5  (curve  2  in  Fig.  1),  and  the  radius  of 
orbit  is  r  =  /?th/^o  multivalued  solution  the  examination  of 

Eqs.  (1.13-1.16),  linearized  in  the  close  vicinity  of  the  steady-state  solutions  (1.17-1.19), 
shows  that  only  those  states  are  stable  which  satisfy  the  energy  criterion 

d(0^)/d{p^)  >  0  (1.21) 

(solid  branches  of  the  curves  in  Fig.  1);  otherwise,  they  become  unstable  (dashed 
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branches  in  Fig.  1).  This  meeis  the  physical  expectations,  and  leads  to  hysteretic 
behavior  of  the  electron  under  conditions  (1.20). 

Let  us  make  some  quantitative  estimates.  A  magnetic  field  of  strength 
Ho  =  100  kG  produces  the  cyclotron  frequency  Uq  =  1.7  x  10*^  sec"*  (Xq  ~  1.07  mm). 
Then  in  the  ultimate  case  of  radiation  losses,  the  resonance  width  is  q  10”**,  which 
yields  the  threshold  field  amplitude  to  be  as  small  as  Etj,  =  1.7  x  10~*  V/cm,  and  the 
kinetic  energy  as  small  as  fig/2  si  1.2  x  10"**.  This  is,  in  fact,  even  near  a"*  times 
smaller  than  a  quantum  limit  of  the  energy  of  excitation  which  is  2!kJo/mc*  (here 
Q  =  e^/lic  =  1/137  is  the  fine-structure  constant).  Therefore,  in  the  close  vicinity  of 
the  threshold  (1.20),  only  the  quantum  approach*^  can  give  an  adequate  description  of 
the  phenomenon,  whereas  for  sufllciently  strong  driving  field  {fi  »  the  classical 
results  (and,  in  particular,  hysteretic  jumps)  remain  valid. 

Very  recently,  the  hysteretic  cyclotron  resonance  of  a  free  electron  in  vacuum  has 
been  first  experimentally  observed  by  G.  Gabrielse  et.  al®  using  a  single  electron  kept 
continuously  in  a  Penning  trap  for  more  than  10  months.  The  hysteresis  was  so  pro- 
noiinced  that  it  provided  the  best  signal  to  noise  ratio  ever  observed  with  a  single  parti¬ 
cle  in  a  trap.  The  cyclotron  energy  was  measured  via  a  second  manifestation  of  the 
relativi.stic  mass  increase.  The  electron  was  weakly  confined  in  a  Penning  trap  and 
oscillated  in  a  direction  which  is  orthogonal  to  the  cyclotron  motion  with  a  frequency 
that  was  measureably  shifted  in  proportion  to  the  electron's  kinetic  energy.  The  typi¬ 
cal  hysteresis  observed  in  the  experiment’^  is  depicted  at  Fig.  2. 

2.  Bistable  cyclotron  resonance  in  serniconductors 

Here  we  shall  consider  a  bistable  (hysteretic)  cyclotron  resonance  in  semiconduc¬ 
tors^.  This  elTect  is  feasible  due  to  the  nonparabolicity  of  the  semiconductor  conduc¬ 
tion  band  which  causes  a  pseudo-relativistic  behavior*’^  of  the  effective  mass  of  conduc¬ 
tion  electrons  in  the  narrow-gap  semiconductors  such  as  c.g.  InSb. 
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The  main  features  of  this  effect  as  compared  with  the  free-electron  case,  are  as  fol¬ 
lows: 

i)  The  nonlinearity  of  conduction  electrons  in  semiconductors  is  many  orders  of 
magnitude  larger  than  the  relativistic  nonlinearity  of  the  free  electron.  This 
allows  one  to  attain  a  fairly  low  critical  pumping  intensity  (although  still  much 
greater  than  in  vacuum)  for  observation  of  the  effect  even  taking  into  considera¬ 
tion  the  much  faster  relaxation  in  semiconductors. 

ii)  The  effective  mass  of  the  electron  in  some  semiconductors  (e.g.  InSb)  is  very  small 
which  results  in  a  considerable  increase  of  the  cyclotron  frequency  (up  to  70-80 
times)  as  compared  to  the  free  electron  for  the  same  magnetic  field  (or  allows  one 
to  correspondingly  decrease  the  required  magnetic  field).  In  the  case  of  InSb,  this 
allows  one  to  use  the  10.6  fim  radiation  of  a  CO2  laser  with  a  magnetic  field  ~ 
140  kG. 

iii)  In  semiconductors,  the  nonlinear  cyclotron  resonance  in  the  strong  optic  resonant 
field  is  accompanied  by  some  small  electrostatic  potential  between  the  faces  of  the 
semiconductor  layer.  This  effect  is  due  to  the  radiation  force  of  the  driving  EM 
field  and  to  the  redistribution  of  electric  charges  in  semiconductor.  This  potential 
exhibits  a  hysteretic  behavior  as  well.  Therefore,  this  effect  could  be  the  first  pro¬ 
posed  all-optical  nonlinear  phenomenon  which  results  simultaneously  in  optic  as 
well  as  opto-electronic  ’  in.tability.  This  property  provides  for  an  immediate  elec¬ 
tronic  indication  of  the  effect. 

We  shall  demonstrate  the  feasibility  of  this  effect  using  a  classical  model  for  the 
interaction  of  an  optical  wave  with  a  single  electron  in  the  conduction  band  of  a  sem¬ 
iconductor.  The  thin  semiconductor  film  is  immersed  in  the  homogeneous  magnetic 
field  Hq  which  is  perpendicular  to  the  layer.  The  semiconductor  is  also  subject  to  the 
action  of  an  EM  wave  of  amplitude  E  which  propagates  along  the  same  direction  as 
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the  magnetic  field  (axis  z),  the  same  as  in  the  free-electron  case.  However,  now  the 
small  static  electric  field  ^  which  is  directed  along  the  z  axis  and  compensates  the  radi¬ 
ation  force,  is  not  provided  anymore  by  the  external  potential,  but  is  caused  by  the 
background  charge.  This  is  due  to  a  redistribution  of  the  electron  density  between  the 
faces  of  the  semiconductor  film,  i.c.  the  problem  becomes  self-consistent  as  far  as  this 
potential  is  concerned. 

In  narrow-gap  semiconductors  which  can  be  described  by  the  Kane  two-band 
model  ®  with  isotropic  nonparabolic  bands,  the  conduction  band  energy  W  (which  is  an 
analog  of  kinetic  energy)  can  be  written  as 

W(p)  =  v/m„*2v<  +  p2v2  (2.1) 

where  ^  is  the  momentum  of  the  conduction  electron,  m^  is  its  effective  mass  at  the 
bottom  of  the  conduction  band,  v^  =  \/VVQ/2m*  is  some  characteristic  speed,  and  Wq 
is  the  band  gap  (the  energy  W  in  Eq.  (2.1)  b  measured  with  respect  to  the  middle  of 
the  gap).  The  velocity  T  of  the  conduction  electron  is  given  by*®  \^(1?)  =  5W(p)/^. 
By  virtue  of  Eq.  (2.1),  this  yields 

-p  -  m^/y/ 1  -  vVv^  ;  t  =  +pVPo^  ,  (2-2) 

where  p^  =  m^Vj,  =  is  some  characteristic  momentum.  One  can  see  from 

Eqs.  (2.1)  and  (2.2)  that  relations  among  W,  v  and  p  are  completely  relativistic,  with  v^ 
playing  a  role  as  an  “effective  speed  of  light,”  and  W(3/2  as  an  “effective  rest  energy”  of 
the  electron.  For  InSb,  Wq  =  0.24  eV,  m*  =  0.014  m^  (where  m^  is  the  rest  mass  of  an 
electron  in  vacuum),  so  that  v^  ~  1.15  x  10®  cm/sec.  Therefore,  it  is  natural  to  look 
for  applications  of  thb  property  of  conduction  electrons  in  solid  state  analogous  to 
those  of  free  electrons.  For  example,  since  the  relativistic  mass-effect  of  the  free  elec¬ 
tron  in  vacuum  is  a  basic  phenomenon  for  the  cyclotron  maser  (“gyrotron”)  ’®'*®,  it  was 
recently  proposed  to  develop  a  solid-state  cyclotron  maser  by  exploiting  the  pseudo- 


relativistic  behavior  of  conduction  electrons  in  semiconductors. 

From  Eqs.  (2.1)  and  (2.2)  one  readily  gets  a  formula  for  the  effective  mass  of  con¬ 
duction  electron  m*(v); 

m*  =  m„*(  I  -  I  v|  (2.3) 

which  again  has  a  relativistic  form,  see  Elq.  (1.3).  Therefore,  one  can  readily  arrive  at 
the  equation  of  the  electron  motion,  Eq.  (MO),  where  now  ^  is  interpreted  as  ^  =  ^/v^, 
u,’o  =  cHo/mgC,  ko  =  a>o\/c(wJ/c,  and  all  the  rest  of  variables  in  Eq.  (1.8)  are  defined 
through  m^  instead  of  m^.  This  immediately  leads  to  the  truncated  dynamic  equations 
(1.13)-(1.16),  and  therefore  -  to  the  equations  (1.17)-(1.19)  for  the  steady-state  regime 
and  the  threshold  conditions  (1.20).  The  only  new  feature  is  that  now,  since  the  quasi¬ 
static  electric  field  g  (or  variable  /?=eg/m*cko)  is  developed  by  the  redistribution  of 
electrons  inside  the  film,  these  equations  are  to  be  supplemented  by  the  equation  for  the 
field  g.  Strictly  speaking,  the  single-electron  model  cannot  provide  the  required  equa¬ 
tion.  The  complete  theory  of  the  dynamics  of  the  system  can  be  developed  only  in  a 
framework  of  a  kinetic  approach  which  takes  into  consideration  the  distribution  func¬ 
tion  of  electrons  in  the  system.  However,  in  order  to  get  results  for  the  steady-state 
regime  as  well  as  equations  for  the  small  perturbations  of  that  regime,  one  can  use 
these  argiiments.  By  the  analogy  with  other  effects  in  which  the  redistribution  of 
charges  is  caused  by  an  external  force  (e.g.  as  in  Hall  effect,  or  in  the  case  when  an 
external  electrostatic  field  is  applied  to  the  film),  it  is  easy  to  understand  that  the  dis¬ 
tribution  of  electrons  remains  almost  homogeneous  inside  the  film  such  that  no  free 
bulk  charges  occur  while  the  induced  charges  are  localized  at  the  surfaces  of  the  film. 
The  thickness  of  the  charged  surface  layers  is  of  the  order  of  the  Debye's  radius  rp, 
which  is  fairly  small  (e.g.  for  <~5  and  cm"®  where  is  the  concentration 

of  electrons,  one  has  rD~10“®  cm  at  room  temperature).  The  entire  film  then  can  be 

considered  as  a  capacitor  with  capacitance  C  =  — ^  per  unit  area  (where  d  is  its 


thickness)  whose  potential  U  determines  the  field  g;  g  =  U/d  =  q/Cd.  Here,  the  electr¬ 
ical  charge  per  unit  area  q  is  related  with  the  current  J  per  unit  area  by  the  relation¬ 
ship  dq/dt  =  J  =  <tE,  where  <t  is  a  conductivity  of  the  material,  and 
Ej  =  \/7E(v/c)cos0-g. 

From  these  relationships  one  gets  q  =  Cdg  =  fg/47r,  or,  after  taking  time  deriva¬ 
tive,  fg/4ff  =  q  =  J  =  <tE,.  Making  use  of  Eq.  (1.15),  one  finally  gets  the  required 
equation  for  g: 

r,,edg/dt  =  m*«v^>  +  =  e  <  E,  >  (2.4) 

where  is  a  relaxation  time  for  the  redistribution  of  the  charge,  and  angle 

brackets  denote  an  averaging  over  ensemble  of  electrons.  The  steady-state  magnitude 
of  g  (or  p),  however,  is  still  determined  by  Eq.  (1.18),  i.e.  or 

egs  =  nio*ckoTf/9*.  (2.5) 

The  threshold  conditions  (1.20)  with  application  to  the  amplitude  of  the  required  criti¬ 
cal  field  amplitude  Etj,  gives 

I  eE|  2  >  I  eEthI  ^  >  qv/S  (2.6) 

Let  us  make  some  qualitative  estimates  for  the  threshold  intensity  Et*  of  the 
incident  wave.  In  the  case  of  InSb,  m^  =  0.014  m<„  Wq  =  0.24  eV.  The  damping  rate 
7  depends  primarily  on  the  concentration  of  impurities.  We  assume  it  to  be  7  —  10“^ 
for  Xo  ~  10.6  fim,  if  a  CO2  laser  is  used  a.*  a  source  of  the  driving  radiation.  The  mag¬ 
netic  field  required  to  get  a  cyclotron  resonance  in  the  neighborhood  of  X  ~  10.6  fim  is 
Ho  ~  140  kG,  and  Eq.  (2.6)  provides  for  a  critical  field  required  to  observe  hysteresis: 
Eih  ~  300  V/cm.  This  corresponds  to  the  incident  power  of  about  240  W/cm*  which 
can  be  readily  attained. 

Let  us  also  estimate  the  quasi-static  potential  U  =  dg  (d  is  the  thickness  of  the 
layer)  which  accompanies  the  phenomenon.  F^om  Eq.  (2.5),  one  can  see  that  the  field  g 


must  experience  hysteresis  as  long  as  the  velocity  0^  does.  At  the  point  of  higher  exci¬ 
tation  where  ~  /i/7  [see  Eq.  (1.13)},  this  gives 

gin»*  -  eE*a;„\/£/Tfm*c  (2.7) 

In  the  above  instance,  if  the  pumping  intensity  is  slightly  above  the  threshold  Ej*  (for 
instance,  E  =  500  V/cm),  one  gets  g„„  •>'  1  V/cm.  For  a  thickness  of  the  layer  d  ~  O.l 
cm  this  gives  U  =  dg  ~  0.1  V  which  can  be  readily  measured.  The  relaxation  time  for 
the  potential  U  is  of  order  of  =  t/iva. 

The  presence  of  two  different  relaxation  times  (r^  and  7“*  •  in  the  system  can 
cause  even  more  complicated  behavior  which  depends  on  the  ratio  of  these  times  i.e. 
quantity  <Jq.  In  the  limit  — ►  0,  i.e.  if  r^7  lJq  <<  1,  it  can  be  shown  that  two 
states  out  of  the  three  possible  steady-states  (in  the  case  of  a  multi-valued  solution)  are 
stable  with  regard  to  small  perturbations.  These  two  stable  steady-states  must  satisfy 
the  same  condition  (1.21)  as  for  free  electrons,  i.e.  those  are  the  states  which  belong  to 
upper  and  lower  branches  of  hysteresis.  However,  if  r^7  Wq  ~  1,  there  could  be  some 
domains  on  the  upper  branch  (which  are  close  to  the  onset  of  the  reverse  jump  from  the 
upper  branch  to  the  lower  one)  where  these  steady  states  becomes  unstable.  This  may 
result  in  excitation  of  self-oscillations  and  possibly  chaotic  motion.  Using  the  usual  for¬ 
mula  for  conductivity  **  (T  =  eNM  where  N  is  the  concentration  of  conduction  electrons 
and  M  =  eujmg'i  is  their  mobility,  one  gets  the  expression  for  parameter 


where  is  the  plasma  frequency  In  the  above  mentioned  instance  of  InSb  with 
7  ^  10"^  this  yields  7  ~  10'*/N  where  N  is  expressed  in  cm"®.  One  must  note, 
though,  that  a  more  comprehensive  theory  of  stability  (as  well  as  of  the  entire 
phenomenon)  has  to  involve  consideration  of  cooperative  motion  of  electrons,  i.e.  the 
kinetics  of  distribution  functions  rather  than  the  motion  of  a  single  electron. 


.  11  - 


3.  Bistable  noncyclotron  resonance  of  a  slightly  -  relativistic  electron 

In  the  above  sections,  the  bistable  elToot  was  discussed  which  is  based  on  a  rotation 
of  a  slightly  •  relativistic  electron  in  a  magnetic  field  (i.e.  on  a  nonlinear  cyclotron 
effect).  However,  it  is  obvious  that  the  presence  of  a  magnetic  field  is  not  a  necessary 
condition;  the  only  factors  substantial  for  the  bistable  resonant  excitation  to  exist  are 
presence  of  nonlinearity  (in  our  case  -  relativistic  mass-efi’ect)  and  sufficiently  sharp 
resonance.  The  latter  one  can  be  provided  by  any  potential  well.  In  the  simplest  and 
probably  most  characteristical  case  it  is  one-dimensional  parabolic  potential  well  (which 
would  correspond  to  a  conventional  linear  oscillator  has  the  electron  mass  not  varied 
due  to  relativistic  effect).  We  assume  that  the  electron  oscillates  along  the  x  axis,  and 
is  excited  by  the  driving  periodic  electric  field  =  E  cos  wt,  which  is  directed  along 
the  same  axis  (for  the  sake  of  simplicity,  we  shall  neglect  here  the  action  of  the  radia¬ 
tion  force,  since  it  does  not  affect  directly  the  steady-state  regime,  see  Section  1),  and 
that  the  motion  is  only  slightly  relativistic  (i.e.  x^  «  c^).  In  such  a  case,  the  motion 
of  the  electron  is  governed  by  the  equation 

X  +  27u;oX  +  xWo(l  -  sin  ut  (3.1) 

nrirt 


where  u)^  is  an  eigenfrequency  of  a  respective  linear  oscillator  in  the  potential  well.  We 
use  the  same  approximation  as  in  Section  I  and  look  for  the  forced  solution  of  Eq  (3.1) 
in  the  form 


x{t)  w  v„,(t)  (jj  *  sin(wt  +  <l>)  (3.2) 

where  the  maximal  speed  of  electron  v,„  and  its  phase  ^  are  slow  varying  functions  of 
time.  We  introduce  again  a  dimensionless  variables 


0m  = 


A  = 


ur-u„ 


-  eE 


1 


m„c- 


(3.3) 


and  a.ssume  that  all  of  them  are  much  smaller  than  unity.  Then,  in  the  same  fashion  as 


in  Section  1,  one  gets  a  truncated  equation  for  their  dynamics: 

+  A*  (t>/2  (3.4) 

H^O  =  "(^  +  i^m/32  +  n  sin  (3.5) 

which  (with  modified  coefficients  and  /?,  =  0)  reproduce  respectively  Eqs  (1.13)  and 
(1.14).  The  steady-state  regimes  (d/dt  =  0)  follows  immediately: 

/|2  =  4^92)^  +  (A  +  0l/32)^  (3.6) 

which  again  reproduces  Eq  (1.17)  with  modified  coefficients.  By  introduction  ji  =  /i/8, 
and  0g  =  ^m/4,  Eq  (3.5)  is  reduced  exactly  to  Eq  (1.17)  for  ji  and  Therefore,  taking 
into  consideration  Eqs  (1.20)  and  (3.3),  the  threshold  conditions  required  in  order  to 
observe  bistability  of  exitation  of  slightly  -  relativistic  oscillator  (3.1)  are 

|eE|  >  32  k„moc2(7/v/3)V2  ;  w  -  a;„  >  v/3  .  (3.7) 

If  the  damping  rate  7  is  again  due  to  radiation  of  EM  wave  by  the  electron,  the 
required  threshold  intensity  of  the  driving  field  is  of  the  same  order  of  magnitude  as 
that  one  discussed  in  Section  1  for  the  bistable  cyclotron  resonance,  in  vaccuum. 

Conclusion 

In  conclusion,  we  demonstrated  the  feasibility  of  hysteretic  behavior  and  bistability 
of  cyclotron  resonance  of  free  electrons  in  vacuum  and  conduction  electrons  in  semicon- 
dtictors  under  the  action  of  sufficiently  strong  quasiresonant  driving  radiation.  We 
showed  also  that  the  same  effect  must  be  peculiar  for  a  conventional  (i.e.  noncyclotron) 
resonance  of  a  slightly  -  relativistic  electron  in  a  one-dimensional  potential  well.  Future 
research  should  involve  quantum  as  well  as  kinetic  theory  of  the  phenomenon.  Even 
far  from  the  onset  of  the  hysteresis,  the  action  of  the  strong  pumping  should  cause  a 
dramatic  change  in  the  location  and  the  shape  of  the  resonant  curve  of  cyclotron 
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resonance;  in  particular,  the  dependence  on  the  frequency  (or  magnetic  field)  should 
become  drastically  asymmetrical.  This  effect  may  provide  a  new  experimental  method 
to  measure  the  nonparabolicity  of  the  conduction  band  in  semiconductors,  effective 
mass,  nonlinear  relarcation,  etc.  Probably,  the  most  attractive  and  fundamental  feature 
of  all  these  efi'ects  is  that  for  the  first  time  they  provide  an  unique  opportunity  to 
experimentally  (and  theoretically)  study  hysteretic  and  bistable  phenomenon  at  a  quan¬ 
tum  level. 
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Captions  to  figures 


The  plots  of  normalized  kinetic  energy  of  the  electron  /^/2  (a)  vs  normal* 
ized  resonant  detuning  A/7  for  various  intensities  of  incident  EM  field,  and 
(b)  vs  normalized  incident  intensity  for  various  detunings.  Curves: 

(a)  1,  /iV2f  =  1;  /iV2y  =  (2>/3f,  3,  1^(2^  =  3;  (b)  1,  A/7  =  0;  2, 
A/7  =  ->/3;  3,  A/7  =  -3. 

Experimentally  observed^  axial  frequency  shift  (left  vertical  scale)  as  a  func¬ 
tion  of  the  frequency  of  the  driving  force  when  this  frequency  is  swept 
downward  (2a)  and  upward  (2b)  through  resonance.  The  right  vertical  scale 
is  the  kinetic  energy  in  the  cyclotron  motion  in  eV. 
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Many  prior  art  bistable  optical  devices  require  resonant 
optical  cavities  and  are  therefore  limited  in  their  opera¬ 
tion  due  to  the  long  lifetimes  associated  with  their  high- 
flnesse  cavities.  A  bistable  optical  device  that  does  not 
use  a  resonant  cavity  is  disclosed  wherein  a  nonlinear 
medium  whose  index  of  refraction  increases  with  in¬ 
creased  light  intensity  is  arranged  to  have  input  and 
output  faces  into  which  and  out  of  which  a  laser  beam 
having  a  nonuniform  spatial  profile  can  be  propagated. 
A  mirror  having  a  predetermined  area  of  reflectivity  is 
positioned  with  respect  to  the  output  face  of  a  nonlinear 
medium  so  as  to  reflect  only  the  light  energy  that  propa¬ 
gates  in  an  area  at  the  output  face  that  is  approximately 
equal  to  the  area  which  the  beam  presents  at  this  face 
when  the  beam  is  propagating  at  a  critical  power  level, 
that  is,  when  the  beam  is  self-trapped. 
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NONLINEAR  OPTICAL  DEVICE  USING 
SELF.TRAPPING  OF  LIGHT 

UACKCROUND  of  the  INVENTION 

The  study  of  optically  bistable  devices  and  their  ge¬ 
neric  characteristics  has  received  increased  attention  in 
the  scientific  community.  See,  for  example,  the  special 
issue  on  optical  bistability  IEEE  Journal  of  Quantum 
Electronics,  QE-17,  March  1981.  Bistable  optical  devices 
are  often  classified  into  one  of  two  categories.  In  the 
first  category  are  intrinsic  devices  which  are  those  de¬ 
vices  in  which  the  feedback  required  for  bistability  is 
optical.  In  the  second  category  are  hybrid  devices  in 
which  devices  some  form  of  electrical  feedback,  some¬ 
times  in  conjunc  n  with  optical  feedback,  is  used. 
Intrinsic  devices  are  of  particular  interest  because  of 
their  potential  for  ultra-fast  switching. 

Many  of  the  intrinsic  devices  known  in  the  prior  art 
require  resonant  optical  cavities.  As  a  result  of  this 
characteristic,  these  prior  art  intrinsic  devices  require 
that  the  input  light  to  the  device  be  tuned  to  a  special 
frequency.  In  addition,  these  devices  are  frequently 
sluggish  in  their  operation  due  to  the  long  lifetimes 
a.ssociated  with  their  high-finesae  cavities.  It  is  expected 
that  the  operation  of  an  intrinsic  optically  bistable  de¬ 
vice  could  he  improved  if  the  device  did  not  require  a 
resonant  cavity. 

SUMMARY  OF  THE  INVENTION 

The  present  invention  is  based  upon  a  fundamentally 
new  type  of  intrinsic  optical  bistability.  The  operation 
of  devices  using  the  present  invention  is  based  on  self- 
focusing  of  light.  Self-focusing  occurs  when  a  light 
beam  having  a  nonuniform  spatial  profile,  such  as  a 
Gaussian  laser  beam,  progagales  through  a  nonlinear 
medium  having  an  index  of  refraction  that  increases 
with  increased  light  intensity.  When  the  light  intensity 
increases  to  a  critical  power  level.  Per,  the  input  laser 
beam  passes  through  the  medium  with  no  change  in 
spot  size,  and  this  siluation  is  referred  to  as  self-trap- 
ping. 

In  accordance  with  the  present  invention  a  nonlinear 
medium  whose  index  of  refraction  increases  with  in¬ 
creased  light  intensity  is  arranged  to  have  input  and 
output  faces  into  which  and  out  of  which  a  light  beam 
having  a  nonuniform  spatial  profile  can  be  propagated. 
A  mirror  is  combined  with  the  output  face  of  the  nonlin¬ 
ear  medium  so  as  to  reflect  only  the  light  energy  that 
propagates  in  an  area  at  the  output  face  that  is  approxi¬ 
mately  equal  to  the  area  which  the  beam  presents  at  this 
face  when  the  beam  is  propagating  at  the  critical  power 
level,  that  is.  when  the  beam  is  self-trapped.  When  the 
beam  power  has  been  increased  to  the  critical  power 
level,  a  substantial  fraction  of  the  beam  will  be  reOecied 
back  into  the  nonlinear  medium.  The  beam  intensity  can 
then  be  decreased  to  a  lower  power  level  and  the  beam 
will  remain  self-trapped  because  of  the  light  energy  that 
is  reflected  from  the  mirror,  and  the  device  will  there¬ 
fore  exhibit  a  bistable  optical  characteristic. 

In  the  embodiment  which  was  constructed,  a  lens  is 
positioned  behind  the  output  face  of  the  nonlinear  me¬ 
dium  in  order  to  image  the  light  at  the  output  face  of  the 
medium  onto  a  light  absorbent  disk  having  an  aperture 
approximately  equal  in  size  to  the  spot  size  that  the 
beam  would  have  at  the  disk  when  the  beam  is  self- 
trapped.  A  mirror  is  positioned  on  the  side  of  the  opti¬ 
cally  absorbent  disk  that  is  opposite  to  the  propagating 
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beam.  The  mirror  is  selected  to  be  partially  transmitting 
in  order  to  provide  an  output  optical  beam. 

BRIEF  DESCRIPTION  OF  THE  DRAWING 

S  The  invention  will  be  more  readily  understood  after 
reading  the  following  detailed  description  in  conjunc¬ 
tion  with  the  drawing  wherein: 

FIG.  1  is  a  schematic  block  diagram  of  an  apparatus 
constructed  in  accordance  with  the  present  invention, 
19  and 

FIG.  2  is  a  graph  of  optical  output  power  versus  input 
power  for  the  device  shown  in  FIG.  1. 

DETAILED  DESCRIPTION 

The  basic  principles  upon  which  the  new  class  of 
device  operates  can  be  understood  from  a  discussion  of 
the  apparatus  shown  in  FIG.  1.  A  laser  beam  10|  of 
power.  P,a,  having  a  Gaussian  intensity  profile  is  fo¬ 
cused  onto  the  input  face  102  of  a  nonlinear  medium  103 
^  having  an  intensity  dependent  index  of  refraction.  For 
self-focusing  to  be  possible  the  medium  must  be  chosen 
such  that  its  refractive  index  increases  with  increasing 
light  intensity.  A  measure  of  the  strength  of  this  nonlin¬ 
earity  is  called  the  “critical  power,"  Per-  When  Pn=Per, 
’  the  input  laser  beam  passes  through  the  medium  with  no 
change  in  spot  size,  as  shown  by  dashed  lines  104  and 
lOS;  this  situation  is  referred  to  as  “self-trapping."  For 
Pui<Ptr  the  laser  beam  diverges  less  rapidly  than  it 
^  would  in  the  absence  of  the  nonlinearity  and  for 
Pi>P<-,  the  beam  converges.  For  P|<<Pf,  there  is 
essentially  no  self-focusing;  the  input  focal  spot  size  is 
chosen  such  that,  for  these  conditions,  the  laser  beam 
diverges  appreciably  in  passing  through  the  nonlinear 
medium,  as  shown  by  the  solid  lines  106  and  107. 

The  optical  field  at  output  face  101  of  the  nonlinear 
.medium  103  is  imaged  by  a  lens  109  onto  the  partially 
transmitting  mirror  111  which  is  aligned  normal  to  the 
laser  beam.  Immediately  in  front  of  the  mirror  is  an 
40  appropriately-sized  optically  absorbent  disk  110  having 
an  aperture  positioned  to  permit  the  self-trapped  beam 
to  pass  through  to  the  mirror  111.  The  aperture  size  is 
adjusted  to  satisfy  two  criteria.  First,  it  must  be  small 
enough  that  in  the  absence  of  self-focusing  the  fraction 
45  of  the  incident  light  feedback  into  the  nonlinear  medium 
by  the  mirror  is  small.  Secondly,  it  is  large  enough  that 
under  self-trapping  conditions  essentially  all  the  light 
passes  through  the  aperture  and  is  reflected  back  upon 
itself  by  the  mirror.  Tliis  strong  feedback  reinforces  the 
50  self-focusing  in  the  nonlinear  medium  and  it  allows 
self-trapping  to  be  maintained  even  when  the  input 
power  is  reduced  below  Per.  This  is  the  mechanism 
which  gives  rise  to  the  optical  bistability  and  hysteresis. 
The  aperture  and  mirror  could  be  combined  and  depos- 
55  ited  on  the  exit  face  of  the  nonlinear  medium  to  mini¬ 
mize  device  length  and  optical  transit  times.  Clearly 
other  arrangements  can  be  devised  by  those  skilled  in 
the  art  which  lead  to  optical  bistability  using  the  same 
basic  principles. 

60  In  the  embodiment  which  was  constructed,  atomic 
sodium  vapor  was  utilized  as  the  nonlinear  medium  103. 
Siodium  vapor  has  the  virtue  of  having  a  large  nonlinear 
index  for  frequencies  near  its  resonance  transitions.  In 
addition,  steady-state  self-focusing  and  self-trapping 
6S  have  been  observed  in  it  using  cw  dye  lasers.  See  the 
article  entitled  “cw  Self-Focusing  and  Self-Trapping  of 
Light  in  Sodium  Vapor,"  by  J.  E.  Bjorkholm  and  A. 
Ashkin,  Physical  Review  Letters,  Vol.  J2,  pages  129-1 JI, 
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Jan.  28,  1974.  The  sodium  vapor  was  contained  in  a 
20-cm*long  heated  cell  120  constructed  of  pyres.  The 
input  laser  beam  101  was  obtained  from  a  single-mode 
cw  ring  dye  laser  (Spectra-Physics  Model  380A).  The 
transverse  mode  of  the  laser  was  TEMqo  (Gaussian 
mode)  and  its  focal  spot  size  on  the  input  face  102  of  the 
medium  103  was  approximately  80  pm;  the  correspond¬ 
ing  confocal  parameter  of  the  laser  beam  was  about  6.8 
cm  so  that  in  the  absence  of  self-focusing  the  spot  size 


As  svill  be  readily  apparent  to  those  skilled  in  the  art, 
a  device  that  exhibits  the  bistable  characteristic  of  the 
type  shown  in  FIG.  2  can  function  as  a  memory  ele¬ 
ment.  For  example,  the  device  can  be  biased  at  an  input 
power  level  of  Ps  at  a  point  substantially  midway  be¬ 
tween  the  two  step  changes  shown  in  FIG.  2.  At  this 
bias  level  of  Ps,  the  output  power  level  is  dependent  on 
the  previous  history  of  input  power  levels.  If  the  input 
power  level  is  increased  to  a  value  in  excess  of  130  mW 


on  the  output  face  108  of  the  medium  was  about  480  •<>  •*>«>  returned  to  the  bias  level  of  Ps  the  device  will 
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pm.  A  linear  polarizer  and  a  quarter  wave  plate  situated 
between  the  laser  and  the  vapor  cell  was  used  as  an 
isolator;  thus  circularly  polarized  light  was  incident 
onto  the  cell.  Lens  109  was  a  16-cm  focal  length  lens 


operate  at  the  point  designated  as  202  in  FIG.  2  and  an 
output  power  level  of  about  1.9  mW  will  be  present.  If 
on  the  other  hand  the  input  power  level  is  decreased  to 
a  value  below  9S  mW  and  then  returned  to  the  bias  level 


positioned  to  image  the  optical  field  at  the  exit  face  108  »  of  ‘levice  will  operate  at  a  point  corresponding  to 


of  the  medium,  with  unity  magnification,  onto  Hat  mir¬ 
ror  111  having  a  reflectance  of  94  percent.  A  disk  110 
having  an  aperture  with  a  diameter  of  approximately 
130  pm  was  placed  several  mm  in  front  of  the  mirror. 

Without  feedback  from  the  mirror  strong  self-focus-  ^ 
ing  and  self-trapping  were  readily  observed  with  ap¬ 
proximately  ISO  mW  of  light  tuned  roughly  1  GHz 
above  the  resonant  frequency  of  the  3S|(F=2)-*3P| 
transition  at  3896  A.  The  sodium  density  was  nominally  „ 
2xl0'2  cm~^.  To  observe  bistability  the  mirror  was 
aligned  normal  to  the  laser  beam  and  the  input  light  was 
amplitude  modulated  at  about  30  Hz  with  a  spinning 
transmission  grating  composed  of  closely  spaced  fine 


201  in  FIG.  2  and  an  output  power  level  of  about  0.3 
mW  will  be  present.  In  this  way,  the  device  remembers 
what  the  previous  input  power  level  has  been  and  can 
therefore  function  as  an  optica)  memory. 

20  What  has  been  describ^  hereinabove  is  an  illustra¬ 
tive  embodiment  of  the  present  invention.  Numerous 
departures  may  be  made  by  those  skilled  in  the  art  with¬ 
out  tteparting  from  the  spiril  and  scope  of  the  present 
invention.  For  example,  if  the  medium  103  has  a  sufli- 

23  cientiy  strong  nonlinearity,  the  entire  device  may  be 
scalable  to  short  lengths  giving  the  potential  for  fast 
response  times.  By  focusing  the  input  laser  beam  in  such 
a  device  to  a  spot  size  on  the  order  of  k,  the  wavelength 
of  the  light,  device  lengths  on  the  order  of  SOX  with 
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wir«.  The  input  power  (P„)  of  beam  lOI  wd  the^wer  corresponding  transit  limes  for  the  light  of  about  0.2 


of  beam  112  passing  through  the  mirror  lit  (Pa,,)  were  should  be  possible, 

monitored  as  functions  of  time  using  photodiodes  and  a  What  is  claimed  is- 

dual-channel  digit^  oscilloscope  (Ni^et  Explorer  II);  I  ^  nonlinear  bistable  optical  device  for  use  with  an 
the  signals  were  also  displayed  and  recorded  as  P,,/  vs.  jjgj,,  having  a  nonuniforro  sfMtia)  profile 

33  comprising  a  nonlinear  medium  (103)  having  an  input 
FIG.  2  presents  data  obtained  whCT  the  parameters  of  fjce  (102)  and  an  output  face  (100),  said  nonlinear  me- 

the  device  were  adjusted  to  exhibit  bistable  behavior;  dium  having  an  index  of  refraction  that  increases  with 

the  laser  was  tuned  roughly  1  GHz  above  the  increasing  light  intensity,  and  means  (109, 110  and  111) 

3S|(F=2)— ‘IPj  transition.  As  the  input  power  was  including  a  reflective  central  aperture  positioned  at  the 

increased  upward,  switching  from  a  low-transmission  40  output  face  of  said  nonlinear  m^ium  for  reflecting  light 
sute  to  a  high-transmission  sute  occurred  at  an  input  tack  into  the  medium  only  from  a  limited  central  area  of 

power  of  about  130  mW.  At  the  switching  point,  the  output  face. 

transmission  abruptly  increased  by  a  factor  of  4.1;  the  g.  a  nonlinear  bisuble  optica)  device  as  defined  in 

rise-time  of  the  switch  was  about  20  ^ssec.  With  input  claim  1  wherein  said  limited  centra)  area  has  a  diameter 

power  decreasing,  a  downward  switch  of  the  same  45  approximately  equal  to  the  diameter  of  an  optical  beam 
speed  occurred  at  about  93  mW  of  input  power.  Optical  propagating  within  said  nonlinear  medium  at  the  self¬ 
pumping  of  sodium  undoubtedly  plays  a  role  in  the  trapping  power  level. 

nonlinearity  of  the  medium  and  the  observed  switching  3.  A  nonlinear  bistable  optical  device  as  defined  in 

time  may  be  more  characteristic  of  optical  pumping  claim  1  wherein  the  means  for  reflecting  light  comprises 

than  of  the  switching  proceu  itself.  The  dashed  lines  30  an  optically  absorbent  disk  (110)  having  an  aperture 
show  the  calculated  low  and  high  transmission  limits,  approximately  equal  to  the  diameter  of  an  optical  beam 

based  on  the  measurement  of  43  mW  incident  onto  the  propagating  within  said  nonlinear  medium  at  the  self- 

0.1 3  mm  aperture  for  an  input  power  of  140  mW  and  no  trapping  power  level,  and  a  partially  transmitting  mir- 

feedback.  The  high  transmission  limit  assumes  100  per-  ror  (111)  positioned  adjacent  to  said  disk  with  the  disk 

cent  transmission  by  the  aperture.  The  upward  switch  33  between  said  medium  and  said  mirror. 


occurs  at  a  power  level  roughly  equal  to  Ptr-  Visual 
observation  of  the  resonance  fluorescence  induced  by 
the  laser  beam  showed  dramatically  that,  at  twitching, 
the  laser  beam  abruptly  changes  from  diverging  propa¬ 
gation  to  what  appeared  to  be  self-trapping.  Under  60 
different  conditions  several  switching  levels  and  several 
hysteresis  loops  were  obtained,  perhaps  corresponding 
to  oscillations  of  the  spot  size  u  the  laser  beam  propa¬ 
gated  through  the  self-focusing  medium. 
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4.  A  nonlinear  bistable  optical  device  as  defined  in 
claim  3  wherein  said  means  for  reflecting  light  further 
includes  a  lens  (189)  positioned  between  said  disk  and 
said  medium. 

9.  A  nonlinear  bistable  optical  device  as  defined  in 
claim  4  wherein  said  nonlinear  medium  is  an  atomic 
sodium  vapor  contained  in  a  heated  cell  (120)  con¬ 
structed  of  pyrex. 
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Extreme-Ultraviolet  and  X-ray  Emission 
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School  of  Electrical  Engineering 
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Weet  Lafayette,  IN  41601 


Introduction 


Fast-moving  electrons  emit  electromagnetic  waves  when  moving  from  one 
medium  into  another  with  a  different  dielectric  constant  [1-3].  This  is  known 
as  transition  radiation  and  was  predicted  by  Ginzburg  and  Frank  |l].  In  a  spa¬ 
tially  periodic  medium  the  waves  emitted  at  different  interfaces  interfere  so 
that  a  resonant  emission  is  obtained  when  the  following  condition  is  satisfied 
(2-4]:  _ 

VT  costf  =  c/v  ~  nX/P  (1) 


jLis 

^  ru4. 

661 


U->:£ 


vCviv 


Abstract 

High-energy  electrons  emit  resonant  electromagnetic  radiation  when  pass¬ 
ing  through  a  spatially  periodic  medium.  It  is  conventionally  assumed  that 
ultra-relativistic  electron  beams  are  required  to  obtain  significant  emission.  We 
demonstrate  theoretically  the  feasibility  of  exploiting  solid-state  superlattices 
with  short  periods  to  obtain  both  spontaneous  and  stimulated  emission  in  the 
far-ultraviolet  and  soft  X-ray  range  using  non-relativistic  beams. 


where  X  is  the  wavelength  of  radiation,  t  is  the  period  of  the  spatially  varying 
dielectric  constant  e  of  medium  (it  is  conventionally  assumed  that  the  varia¬ 
tions  are  very  small),  v  is  velocity  of  the  electrons  (assumed  normal  to  the 
interfaces),  $  is  the  angle  between  the  direction  of  wave  propagation  and  elec¬ 
tron  motion,  n  is  an  integer,  and  T  is  a  “mean”  c.  This  condition  is  readily 
derived  by  requiring  that  the  waves  emitted  at  different  interfaces  interfere 
constructively  at  a  distant  point.  Usually  the  period  P  >>  X,  so  that  ultra- 
relativistic  beams  (v/c  ot  1)  are  required  in  order  to  satisfy  Eq(l)  for  real  0. 
Recently  the  possibility  of  stimulated  resonance  radiation  of  ultra- relativistic 
(~  50  Gev)  electrons  traveling  through  a  stack  of  metal  foils  was  considered  [4j, 
with  P  ~  7  cm. 


In  this  paper,  we  demonstrate  the  feasibility  of  using  non-relativistic  elec¬ 
tron  beams  in  order  to  attain  both  spontaneous  and  stimulated  emission  in  the 
ultraviolet  and  X-ray  range  using  solid-state  superlattices  with  P  ~  100  A  so 
that  P/nX  ~  1.  We  show  that  the  wavelength  of  resonant  radiation  and  the 
required  energy  of  electrons  are  determined  by  the  parameter  Q  =  nXp/P, 
where  Xp  is  a  “mean”  plasma  wavelength  of  the  medium.  If  nXp  <<  P  (i.e.  Q 
<<  1),  as  is  assumed  in  all  previous  work  [2-4],  then  the  wavelength  of  the 
resonant  radiation  has  an  order  of  magnitude  of  X  X.Q,  and  the  kinetic 
(dimeasionless)  energy  of  the  electrons  eU/mc^  must  exceed  the  critical  amount 
~l/vQ*-tf*  >>  1  (tf  <  Q)  which  constitutes  the  use  of  ultra-relativistic 
beams.  On  the  contrary,  if  the  period  of  the  spatially  periodic  medium  is 
chosen  so  that  Q  >>1  (i.e.  period  P  is  much  shorter  than  the  plasma 
wavelength,  e.g.  P  ~  50*200  A  ),  the  wavelength  of  resonant  radiation  becomes 
of  the  order  of  ~  Xp/Q  =  P/n,  and  the  criticsl  kinetic  energy  of  the  beam  turns 


out  to  be  extremely  small:  eU/mc^  ~  1/2Q^,  which  may  be  less  than  a  few  kilo¬ 
volts  even  for  very  short  wavelength  radiation.  The  advantages  of  the  pro¬ 
posed  method  are;  (1)  the  frequency  of  radiation  can  be  easily  tuned  in  a  very 
wide  range  by  simply  varying  accelerating  potential  of  beam  (which  is  very 
hard  to  do  with  ultra-relativistic  beams)  m  the  range  of  the  possible  angles  of 
the  wave  propagation  is  almost  unlimited,  and  (3)  the  cost  of  equipment  and 
energy  required  for  experiments  with  non-relativistic  beams  is  insignificant 
compared  to  large  accelerating  machines.  This  last  consideration  is  perhaps  the 
most  important. 

The  main  requirements  for  non-relativistic  short-wavelength  resonant  radi¬ 
ation  is  a  periodic  medium  with  a  very  short  spatial  period.  Fortunately,  the 
development  of  molecular  beam  epitaxy  (MBE)  and  other  techniques  in  recent 
years  has  made  it  possible  to  grow  very  tnin  films  (~  100  A  and  less)  with  pre¬ 
cise  boundaries.  Periodic  structures  composed  of  thin  films  of  different  materi¬ 
als,  in  particular  superlattices,  have  also  been  fabricated  (5).  Using  these  struc¬ 
tures,  non-relativistic  electron  beams  with  energies  20-200  KeV  can  be  used  to 
generate  radiation  of  wavelength  100-200  A  and  less.  The  concept  of  EM  radia¬ 
tion  of  electron  beam  in  spatially  periodic  structures  in  general  is  known  since 
microwave  traveling  wave  tube  amplifiers  (with  nonrelativistic  beams)  and  their 
recent  optical  modification  -  free-electron  lasers  (with  relativistic  beams).  An 
important  feature  of  transition  radiation  discussed  in  this  paper  is  that  the 
electron  beam  travels  through  the  material  structure  rather  than  in  vacuum 
above  the  structure  las  in  traveling  wave  tubes),  or  through  a  spatially  modu¬ 
lated  magnetic  field  (as  in  many  free-electron  lasers).  This  makes  it  possible  to 
use  very  short  spatial  periods.  The  problem  is  the  energy  transfer  from  the 
electron  beam  to  the  material  structure  causing  heating  and  possible  damages. 
This  problem  will  be  briefly  addressed  at  the  end  of  this  paper. 

In  this  paper,  we  show  that  resonant  spontaneous  emission  with  a  total 
power  of  0.1  mW  (around  a  wavelength  ~200  A)  can  be  obtained  with  a  75 
KeV  electron  beam  carrying  a  current  of  only  1  mA.  The  spontaneous  emission 
can  be  used  as  a  narrowband  source  by  selecting  the  radiation  in  a  narrow 
angular  range.  To  get  stimulated  emission  with  a  aain  of  5%  per  pass,  how¬ 
ever,  requires  a  significantly  larger  current  (5xl0‘^A/cm^);  to  avoid  sample 
burnout  it  will  be  necessary  to  use  pulsed  operation  witn  pulse  lengths  M).l  ps, 
which  is  a  difficult  task.  However,  it  should  be  noted  that  coherent  sources  are 
not  available  at  these  short  wavelengths  and  getting  significant  stimulated 
emission  is  a  difficult  problem  in  general. 

Spontaneous  Transition  Emission 

We  will  fir^  obtain^e  resoi^t  wavelength  of  radiation  from  Eq(l)  not¬ 
ing  that  [3)  Jvt  =  Ji\/q  +  where  0|2  and  C|2  are  respectively  the 

thicknesses  and  dielectric  constants  of  alternating  layers  forming  the  superlat¬ 
tice,  9  =01+^2  spatial  period.  For  short  wavelengths  (X^  <<  X|‘2)i 

C]2  =  X|2  are  the  plasma  wavelengths  of  the  two  materi’ab 

forming  the  superlattice,  (Xf2  ~  4]rmc^/e^N|2  where  N*  is  the  density  of  elec¬ 
trons).  Thus,  the  mean  dielectric  constant  may  be  written  as 
e  tat  I  -  XVX*,  where  X"*  =  P“‘(^iXf*+j2^2*)J-  Substituting  this  into  (1)  and 
solving  it  for  X,  one  gets  the  resonant  wavelengths: 

"  =  (2) 
where  Q  =  nXp/f  ;  ^  =  v/c;  7  =  and  =  (H-(Q*-sin^tf)“*j*^*  ** 


critical  energy 


for  the  excitation  of  resonant  radiation.  For  Q  <  <  1, 


X  a<  Xp(Q  ±  Here  Tfer  “  >>  1»  such  that  only  an 

ultra-relativistic  beam  can  excite  radiation.  On  the  other  hand,  when  Q  >> 
1,  the  critical  kinetic  energy  turns  out  to  be  extremely  low, 
(eU/mc*)jj  =  %t~l  «  1/2Q*,  which  is  less  than  10  KeV  for  all  conventional 
materials  if  t  ~  100  A.  For  sufficiently  higher  (but  still  non-relativistic)  ener¬ 
gies  elJ,  Eq(2)  eives  simply 

X  -(“cosfl)  -(\/mcV2eU-cosfi);  eU  <<  me*  =  0.51  MeV.  (3) 
D  jj  a 

For  instance,  if  f=100  A,  n  =  1,  eU  =  75  KeV,  and  tf  =  45“,  one  has 
X  =  113.6  A;  for  n  =  10,  X  =  11.3  A. 

The  resonant  radiation  (i.e.  spontaneous  emission  in  the  system)  can  pro¬ 
vide  a  narrowband  source  of  radiation.  A  single  electron  traversing  multilayer 
structure  with  N  layers  radiates  energy  I  in  a  solid  angle  dO  in  the  frequency 
interval  between  uf  and  cj-hdtj,  given  by  [1-4] 

d^I/dwdn  =  (d\^wdn)  •  4sin2(^P  ^/t  )sin*(^/2) /sin*^  (4) 

where  (  =  (-^-\/i'cosfl)irf/X,  and  d^Io/dfldw  is  a  radiation  produced  by  a  single 

P 

interface.  According  to  Ginsburg-F rank  theory  (1-3].  for  non-relativistic  elec¬ 
trons  <<  1)  and  small  variations  of  c  (j  Afj  =  ^  distri¬ 

bution  of  single-interface  radiation  is  given  by 

d^Io/dfldw  =  e*)9*(Af(w)psin*tf/4Ji^c  (5) 

If  the  number  N  of  layers  is  sufficiently  large,  (N  >>  |c/Ac| ),  Eq(4)  provides 
for  very  narrow  spectral  peaks  of  radiation  for  each  particular  angle  9  (with 
central  wavelength  determined  by  Eq(3)),  which  also  implies  that  any  frequency 
is  radiated  in  a  very  narrow  intervals  of  angles.  Noting  that 
Ac  =  X*(Xf*-X2*),  and  integrating  Eq(4)  over  w  and  fl  (with  dfl  =  2xsin0d0), 
one  gets  the  total  radiation  in  each  order  n 

I  ry  16e*Lj 2(X,-2-X2  2)2sin2(n;r(l ,/J  )/3^nV  ,  (6) 

(where  L  =  N(l/2  is  the  total  thickness  of  the  structure)  with  the  wavelengths 
of  radiation  being  in  the  range  ~(~~1)  <  X  <  “("r  +  l).  The  total  energy 
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of  radiation  increases  as  speed  decreases.  In  order  to  calculate  the  power  of 
resonant  radiation  emitted  by  an  electron  beam  with  an  electrical  current  J, 
one  must  multiply  Eqs(4)-(6)  by  J/e.  If  P  =  100  A,  Pi  =  P2  =  P/2,  L  =1  ftm, 
eU  =  75KeV,  J  =  1mA,  and  Xj  ai  400  A  (e.g.  Zn,  Cu,  Ag  or  Au),  X2  m  800  A 
(e.g.  Si  or  Ge,  see  (6j),  the  system  can  provide  a  radiation  of  first  harmonic  (n 
=  1)  with  a  total  power  0.1  mW  and  a  mean  wavelength  ~  200  A. 


Stimulated  Emiulon  (Amplifleatlou) 

We  will  derive  now  an  amplification  caused  by  the  stimulated  emission. 
This  effect  may  be  viewed  in  the  following  way.  An  EM  wave  having  a  wave 
vector  component  k,  =  k^  cos9  along  the  axis  z  (which  coincides  with  the  elec¬ 
tron  trajectory)  produces  the  higher  order  spatial  harmonics  with 
k,  =  k,  ±  2]ra/P  which  is  due  to  the  periodicity  of  medium.  The  phase  velo¬ 
cities  oL  these  harmonics  along  the  axis  x  are,  therefore, 
~  c/v^cosP^xn/PkJ.  If  the  resonant  condition  (1)  for  X  is  fulfilled,  one 
of  these  phase  velocities  coincides  with  the  speed  of  the  electron  that  results  in 
an  exchange  of  energy  between  the  EM  wave  and  the  electron.  For  some 
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frequencies  in  the  neighborhood  of  resonance,  the  electron  loses  energy  to  the 
EM  wave;  this  results  in  a  coherent  gain  of  the  wave,  or  stimulated  emission. 

Essentially,  this  resembles  a  common  mechanism  of  amplification  for  many 
kinds  of  microwave  devices  based  on  the  interaction  of  electrons  with  "slow” 
EM  wave.  The  important  point  is  to  find  the  intensities  of  the  resonant  spatial 
harmonics  of  the  field.  In  all  the  previous  work  on  resonant  radiation  (2>4]  it  is 
assumed  that  X  <  <  f  (  which  is  always  valid  in  the  ultra^relativistic  case,  see 
the  introductory  section).  This  allows  one  to  use  the  WKB  approximation. 
This  approximation  is  not  valid  in  our  case  since  X  may  be  of  order  or  longer 
that  9.  Instead  we  will  find  a  solution  of  the  exact  wave  equation  (with 
periodic  parameters)  based  on  the  assumption  of  smallness  of  variations  of  sus¬ 
ceptibility  (i.e.,  I  Ac/c  I  <  <  1,  which  is  always  true  for  short  wavelengths);  no 
assumption  is  made  regarding  the  ratio  X/d.  Furthermore,  the  spatial  variation 
of  e(z)  is  usually  approximated  by  a  cosine  function  (2-4).  In  this  approach,  the 
relative  amplitude  of  n^  harmonic  of  the  EM  field  is  p.  ~  (Ae/f)'*,  so  that 
for  small  Ae  p„  is  negligible  for  all  but  the  smallest  n  (=  In  our  approach, 
we  can  treat  any  arbitrary  function  of  ((z),  in  particular  the  true  rectangular 
function.  We  show  that  p„  falls  off  algebraically  like  the  Fourier  coefficients  of 
c(z).  Significant  radiation  is  expected  even  for  large  n  provided  the  interfaces 
are  sharp  enough.  In  this  paper,  we  approach  the  problem  using  a  single¬ 
particle  picture  which  provides  direct  insight  into  the  mechanism  of  the 
electron-EM  wave  interaction.  The  problem  can  also  be  treated  using  either 
the  Boltzmann  equation  [4]  or  a  quantum  mechanical  formalism;  we  plan  to 
address  these  aspects  in  a  subsequent  publication  (7). 

We  consider  the  exact  Maxwell  equation  for  the  EM  field  (9)  with  £(z) 
being  an  arbitrary  periodic  function  in  z.  We  assume  a  plane  wave;  it  can  be 
shown  that  only  the  EM  wave  with  its  electric  field  S  polarized  in  the  plane  of 
incidence  (i.e.  plane  x,z)  may  be  amplified  by  the  beam  [8].  By  virtue  of  the 
Floquet’s  theorem  for  wave  equations  with  periodic  coefficients  [9],  any  com¬ 
ponent  of  the  EM  field  can  be  written  as  a  sum  of  spatial  harmonics; 

u  =  u,exp(jir*P-ju/t)(l  +  X;p„exp(2jnTz/d  +  j^J);  (7) 

n^O 


where  k*T  =  kgXsind+koZcosd,  and  p„  is  the  amplitude  of  the  n^  spatial  har¬ 
monics.  We  make  the  conventional  assumption  that  there  is  no  retrorefiection, 
which  is  valid  if  N(  At/t]  *  <<  1,  and  j  Ae/(J  <<  costf.  This  assumption  is 
strictly  true  in  the  vicinity  of  fi  =  45"  (see  [8]).  Substituting  the  EM  field  in 
the  jprm  (7)  into  the  Maxwell  equations,  collecting  the  terms  with 
exp(jk’T+2jn]rz/(l)  for  each  particular  n  and  retaining  only  terms  that  are  first 
order  in  ai],|<<cos^),  where  the  a^'s  are  the  Fourie  coefficients  of  e(z); 

£(z)  =  £  +  5]  a„cos(2nirz/f  +V’ii)i  8®^  ^l*®  amplitudes  of  the  spatial 
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harmonics  of  nonvanishing  con^imnts  of  electric  and  magnetic  fields 
(Ej|,Hy,E,):  E,  =  E„costf;  H,  =  E„V£;  E,  =  -E„sintf,  and 

^  O 


'0 

1  +qsin*tf/co8tf 

^i/2 

1+qcostf 

q(2cos9-t-q) 

Pt 

a 

2 

l-qcostf-q 

(8) 


where  E^  is  the  ampli^de  of  the  principal  harmonic  of  total  electric  field,  and 
q  =  2xD/tk„  =  \d/9vT.  Further  calculations  are  based  on  the  conventional 


^  •. 
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model  of  energy  exchange  between  the  EM  field  and  an  electron  which  is  used, 
e.g.  in  the  theory  of  free-electron  lasers  (see  e.g.  [10]).  From  the  Lorentz  equa¬ 
tion  me  d(^)/dt  =  e(S+|^H|).  one  gets  the  equation  for  the  energy 
(  ~  Tfmc^  of  electron  -  _  . 

di/dt  =  ec(^)  =  ec(^B+lBA]5+^AB);  (9) 

where  £  =  E![F(t),t|  is  the  field  at  the  instantaneous  location  of  the  electron,  ^ 
and  S  are  unperturbed  v^tors,  A^S  is  a  small  perturbation  of  electron  velocity 
due  to  interaction,  and  AE  is  a  small  perturbation  of  the  field  seen  by  the  elec¬ 
tron  due  its  spatial  displacement  in  respect  to  the  unperturbed  trajectory,  i.e. 

aff  * 

AE  =  —  Az(t);Az  =  c  /  A/d,dt  (10) 

dz  5 


For  the  assumed  polarization  of  the  field,  it  follows  from  the  Lorentz  equation 
that 

t  t 


=  ^  /  E.dt 

7***®  0  TT me  0 


(11) 


In  Eqs  (9)-(ll)  one  has  to  take  into  account  only  that  particular  n^  component 
of  the  wave  which  b  “resonant”  to  the  speed  of  electron,  i.e.  that  one  with 
I  (~n]r|  <  <  1.  After  substituting  z  =  and  the  amplitudes  (8)  of  the 
proper  resonant  harmonic  of  E,,Hy,  and  E,,  into  (9)-(ll),  integrating  over  the 
temporal  interval  |0,r  =  L//9cl,  wnere  r  is  a  time  for  an  electron  to  pass 
through  the  superlattice,  one  has  to  average  the  result  over  all  the  possible 
phases  of  the  relevant  field  harmonic.  We  denote  this  operation  by  angle 
brackets.  Note  that  the  term  <^E>  in  (9)  vanishes,  i.e.  the  stimulated  emis¬ 
sion  is  only  due  to  changes  in  the  electron  motion  caused  by  the  field.  Finally, 
one  gets  the  total  averaged  change  of  the  electron  energy  per  pass; 

e%^L^  1  1/ 

<A#„>  =  JT  "V/  *  ■; — :r  cose{p,,  cose-^y/zp  )  —  x 


m)9®c^TfX  (i/„T)®  *11 


u 


(l-cosi/nf) +7  Vi*  sin*^-2 +2cosj/„r+t'„7sini/„r)|;  (12) 

8 

where  i/„  =  w[/jVii(costf+q)-lj  =  fic(nir-0/9  is  a  resonant  factor.  For  some 
the  change  of  energy  <AI>  becomes  negative  which  constitutes  the  gain  of 
EM  field,  <A^em>  “  -<AI>.  In  the  non-relativistic  case,  the  main  contri¬ 
bution  to  the  change  of  energy  is  due  to  z-components  of  A^  and  AE  i.e.,  in 

(12)  (p*,  p„py)  <<  p*  M  (a„/2)*(l-/3cosfl)*.  Replacing  the  term 

(-2+2cosi/T+i/7sini/T)/i/*r*  by  its  negative  extremum  -4/jt^  (i/r  »  jt),  one  gets 
the  maximal  EM-wave  gain  per  electron  per  pass: 

<A^EM>  "  sin*tf(l//?-’C08^)V»*nic*/SX  (13) 


In  order  to  obtain  an  amplification  F  per  pass  in  the  system  bombarded  by  an 
electron  beam  with  the  density  of  electric  current  i(A/cm^),  one  has  to  multiply 
<A^Em>  by  i/e  and  divide  by  the  energy  flux  of  incident  EM  wave  per  unity 
ares  of  the  interface  E^costf/2R,  where  R  =  3770  is  the  vacuum  impedance. 
One  has  also  to  t^e  into  account  that  for  rectangular  form  of 
€(z),  an/2  =  (Ac/nir)  sin(n]r(i/()  with  Ac  =  X^jXi^-Xf*).  Bearing  in  mind  a 
resonant  condition  (3),  one  finally  gets  the  maximal  EM  wave  amplification  per 
pass; 

r  =  8pieRL*sin*(jrafi/J)sin*^/mc*f  jr^cos^,  (14) 

where 


H  =^*f<(X,-2-X22)2{l//?-COstf)5Q-5  =X*r‘(Xf2-X2*)2(cOStf+nX/f). 

If  X  =  140A,  t  =  lOOA,  9i  =  t2  =  SOA,  n  =  l,  L  =  l/im,  X,  400  A, 
X^  a«  800  A,  =  45  * ,  and  i  =  5x10*®  A/cm^  f4j  (i.g.  beam  of  2  /im  diameter 
with  a  current  ~  1.5x10^  A),  one  gets  an  ampiincation  P  at  5%  per  pass.  The 
required  speed  of  electrons  is  as  0.474  which  corresponds  to  energy  eU  =  69 
KeV.  For  larger  X  the  amplification  increases  drastically.  With  the  mirrors 
situated  outside  the  superlattice  to  form  a  Fabry-Perot  resonator  to  provide 
feedback,  the  system  becomes  a  short-wave  laser  which  may  transform 
significant  portion  of  energy  of  electron  beam  into  coherent  radiation.  It  is 
obvious  that  the  amplifiers  and  lasers  based  on  the  proposed  principle  should 
work  in  the  short  pulse  regime  of  operation,  with  the  duration  of  current  pulse 
being  determined  by  the  heating,  ionization,  diffusion  of  absorbed  electrons, 
etc.  As  a  rough  approximation,  the  per  atom  heating  rate  caused  by  the 
energy  losses  of  the  electron  beam  [11],  is 

(n*/N*)df/dt  =  4iriZ(mv^)'‘ln('7^mv*/e*‘*^)  (15) 

where  n‘  =  i/ev  is  the  electron  beam  density,  N*  is  the  atomic  density  of  the 
material,  Z  is  the  atomic  number.  For  the  parameters  mentioned  above,  the 
duration  of  the  current  pulse  must  be  shorter  than  ~10~*^  sec  in  order  for  the 
energy  transfer  per  atom  to  be  of  order  of  ~1  eV  or  less.  One  may  note  though 
that  in  the  case  of  ultra-relativistic  beams  |4|  with  energy  '^50  GeV  for  the 
same  current,  the  losses  1 15)  are  even  greater,  such  that  one  needs  even  shorter 
pulses.  It  will  probably  oe  hard  to  obtain  such  short  and  powerful  pulses,  but 
it  may  prove  worthwhile.  However,  the  first  step  is  to  attain  spontaneous 
resonant  emission  as  described  in  the  beginning  of  this  letter.  The  spontaneous 
emission  can  be  used  as  a  verv  narrowband  source  of  radiation  by  selecting  nar¬ 
row  range  of  angles  [Eq.  (4)|.  It  should  also  be  noted  that  the  spontaneous 
radiaton  intensity  depends  on  the  total  current  in  the  electron  beam  unlike  the 
stimulated  emission  gain  which  depends  on  the  current  density.  Since  there  is 
no  constraint  on  the  current  density,  the  conditions  to  obtan  spontaneous  emis¬ 
sion  are  much  more  relaxed.  According  to  Eqs.  (4-6),  the  electrical  current 
reauired  to  observe  spontaneous  emission  with  the  same  spotsize  of  the  beam  is 
10®-10®  less  than  that  required  for  stimulated  emission. 

Conclusion 

In  conclusion,  we  have  demonstrated  the  feasibility  of  generating  far- 
ultraviolet  and  soft  X-ray  radiation  by  electron  beams  with  relatively  low, 
non-relativist ic  energies,  traversing  the  solid-state  superlattice  composed  of 
very  thin  periodic  layers.  The  use  of  low  energies  is  a  desirable  feature  as  com¬ 
pared  with  ultra-relativistic  beams.  The  proposed  system  can  be  used  as  a  very 
efl[icient  noncoherent  source  of  narrowband  radiation,  and,  under  special  condi¬ 
tions,  as  an  amplifier  and  laser. 

We  thankfully  appreciate  very  useful  discussions  with  G.  Ascarelli,  R. 
Gunshor,  as  well  as  with  P.  Kelley,  A.  Calawa,  N.  Economou,  and  all  partici¬ 
pants  at  the  seminar  given  by  A.  E.  K.  at  Lincoln  Lab,  M.I.T. 
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High-energy  electrons  emit  resonant  electromagnetic  radiation  when  passing  through  a  spatially 
periodic  medium.  It  is  conventionally  assumed  that  ultrarelativistic  electron  beams  are  required 
to  obtain  significant  emission.  We  demonstrate  theoretically  the  feasibility  of  exploiting  solid- 
state  superlattices  with  short  spatial  periods  to  obtain  both  spontaneous  and  stimulated  emission 
in  the  extreme-ultraviolet  and  soft  x-ray  range  using  nonrelativistic  beams. 

PACS  numbers:  78.45.  -)-  h,  79.20.Kz.  41.80.Dd 


Fast-moving  electrons  emit  electromagnetic  waves 
when  moving  from  one  medium  into  another  with  a  different 
dielectric  constant.*'^  This  is  known  as  transition  radiation 
and  was  predicted  by  Ginzburg  and  Frank. '  In  a  spatially 
periodic  medium  the  waves  emitted  at  different  interfaces 
interfere  so  that  a  resonant  emission  is  obtained  when  the 
following  condition  is  satisfied'^: 

co&d  =  c/v  —  nA /I,  (1) 

where  A  is  the  wavelength  of  radiation,  /  the  period  of  the 
spatially  varying  dielectric  constant  £  of  medium  (it  is  con¬ 
ventionally  assumed  that  the  variations  are  very  small),  v 
velocity  of  the  electrons  (assumed  normal  to  the  interfaces),  6 
the  angle  between  the  direction  of  wave  propagation  and 
electron  motion,  n  an  integer,  and  ?  a  '‘mean”  e.  Usually  the 
period  />/i,  so  that  ultrarelativistic  beams  (t)/c~l)  are  re¬ 
quired  in  order  to  satisfy  Eq.  (1)  for  real  0.  Recently,  the 
possibility  of  stimulated  resonance  radiation  of  ultrarelati¬ 
vistic  ( ~  50  GeV)  electrons  traveling  through  a  stack  of  met¬ 
al  foils  was  considered,^  with  /~7  cm. 

In  this  letter  we  demonstrate  the  feasibility  of  using 
nonrelativistic  electron  beams  in  order  to  attain  both  sponta¬ 
neous  and  stimulated  emission  in  the  extreme  ultraviolet  and 
x-ray  range  using  solid-state  superlattices  with  100  A  so 
that  /  /nA  —  1 .  We  show  that  the  wavelength  of  resonant  radi¬ 
ation  and  the  required  energy  of  electrons  are  determined  by 
the  parameter  Q  =  nA^/l,  where  A^  is  a  "mean”  plasma 
wavelength  of  the  medium.  If  Q<1,  as  is  actually  assumed  in 
all  previous  work,^“*  then  the  wavelength  of  the  resonant 
radiation  has  an  order  of  magnitude  of  A~A^Q,  and  the  ki¬ 
netic  energy  of  the  electrons  e  U  /mc^  must  exceed  the  critical 
amount  ~ \/\Q^  —  d^>\(6 < Q )  which  constitutes  the  use 
of  ultrarelativistic  beams.  On  the  contrary,  if  the  period  of 
the  spatially  periodic  medium  is  chosen  so  that  Q>I,  the 
wavelength  of  resonant  radiation  becomes  of  the  order 
~Ai,/Q  =  /  /n,  and  the  critical  kinetic  energy  of  the  beam 
turns  out  to  be  extremely  small;  eU /mc^  ~  1/20  ’.  The  ad¬ 
vantages  of  the  proposed  method  are  (1)  the  frequency  of 
radiation  can  be  easily  tuned  in  a  very  wide  range  by  simply 
varying  accelerating  potential  of  beam  (which  is  very  hard  to 
do  with  ultrarelativistic  beams),  (2)  the  range  of  the  angles  of 
the  radiated  emission  is  almost  unlimited,  and  (3)  the  cost  of 
equipment  and  energy  required  for  experiments  with  nonre¬ 
lativistic  beams  is  insignificant  compared  to  large  accelerat¬ 


ing  machines.  This  last  consideration  is  perhaps  the  most 
important. 

Fortunately,  the  development  of  molecular  beam  epi¬ 
taxy  (MBE)  and  other  techniques  in  recent  years  has  made  it 
possible  to  grow  very  thin  films  ( ~  100  A  and  less)  with  pre¬ 
cise  boundaries.  Periodic  structures  composed  of  thin  films 
of  different  materials  (in  particular,  superlattices)  have  also 
been  fabricated.*  Using  these  structures,  nonrelativistic  elec¬ 
tron  beams  with  energies  70-200  keV  can  be  used  to  generate 
radiation  of  wavelength  100-200  A  and  less.  An  important 
feature  of  transition  radiation  is  that  the  electron  beam  tra¬ 
vels  through  the  material  structure  rather  than  in  vacuum 
above  the  structure  (as  in  traveling  wave  tubes''),  or  through  a 
spatially  modulated  magnetic  field  (as  in  many  free-electron 
lasers’).  This  makes  it  possible  to  use  very  short  spatial  per¬ 
iods. 

We  will  first  obtain  the  resonant  wavelength  of  radi¬ 
ation  from  Eq.  ( 1 )  noting  that '  /  v?  =  /|vei  +/:ve;,  where 
l,  i  and  e,  2  are,  respectively,  the  thicknesses  and  dielectric 
constants  of  alternating  layers  forming  the  superlattice; 
l=li+ 1^  is  its  spatial  period.  For  short  wavelengths 
[A  f  2),  c,  2  =  \  —A'/A  f  2-  where  /1, 2  are  the  plasma 
wavelengths  of  the  two  materials  forming  the  superlattice 
{A  f  2  =  ^trmc^/e'N  '2 .  where  A  ■■  is  the  density  of  electrons). 
Thus,  the  mean  dielectric  constant  may  be  written  as 
f~l  —  A  VAp,  where  A  =  l~  '\I^A  +  l^A  ^  ‘I-  Sub¬ 

stituting  this  into  Eq.  (1)  and  solving  it  for  A,  one  gets  the 
resonant  wavelengths: 

A=A  g  r  1  cosg 

"  0"-)- cos’ (9  1)9  0 

X  V(K, -l)-'-(?^-l)-'],  (2) 

where  0  =  v/c,  y  =  (l  — )3')~',  and 
=  (I  +  (0’  —  sin^0)~  is  thecritical  energy  required  for 
the  excitation  of  resonant  radiation.  For  0<1,  Aa^ryi^ 
(Q  ±  ylVcr  ^  Here,  ?'„>!,  such  that  only  an  ultrare¬ 

lativistic  beam  can  excite  radiation.  On  the  other  hand,  when 
0>1,  the  critical  kinetic  energy  turns  out  to  be  extremely 
low,  y„  —  l~l/20^  which  gives  less  than  10  keV  for  all 
conventional  materials  if  /~  100  A.  For  sufficiently  higher 
(but  still  nonrelativistic)  energies  eU,  Eq.  (2)  gives  simply 
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ABSTRACT 

Resonant  transition  radiation,  generated  by  high  energy  electron  beams  traversing 
a  periodic  medium,  has  been  considered  by  many  researchers  as  a  potential  source  of 
both  spontaneous  and  stimulated  emission  at  short  wavelengths.  To  our  knowledge, 
this  problem  has  only  been  treated  cla.ssically.  This  paper  presents  a  quantum  mechan¬ 
ical  theory  that  leads  to  a  unified  description  of  both  spontaneous  and  stimulated  emis¬ 
sion  and  establishes  a  simple  relation  between  them. 


1.  Introduction: 


Electrons  traveling  at  high  speed  emit  electromagnetic  waves  when  they  move 
from  one  medium  to  another  with  a  different  dielectric  constant.  This  is  known  as 
transition  radiation'.  In  a  spatially  periodic  medium  there  is  an  interference  of  the 
waves  emitted  at  different  interfaces  producing  a  resonant  transition  radiation  when  the 
following  condition  is  satisfied  (Fig. 

=  (1) 

where 

X  is  the  wavelength  of  radiation  in  free  space 

9  is  the  period  of  the  spatially  varying  dielectric  constant  E(z)  [the  variations  are 
usually  assumed  to  be  small], 

£  is  the  ‘mean’  dielectric  constant, 

V  is  the  velocity  of  the  electrons  in  the  z-direction  perpendicular  to  the  interfitces, 

&  is  the  angle  between  the  direction  of  wave  propagation  and  the  z*axis, 

and  n  is  an  integer 

Resonant  transition  radiation  has  been  considered  by  a  number  of  workers  as  a 
possible  source  of  short  wavelength  radiation.  Usually  the  period  9  is  much  greater 
than  X  so  that  ultra-relativistic  electron  beams  with  v/c  tst  1  are  required  to  satisfy  Eq. 
1^.  The  possibility  of  using  non-relativistic  electron  beams  with  short  period  solid  state 
superlattices  has  also  been  considered  recently^. 

To  our  knowledge  the  problem  of  resonant  transition  radiation  has  only  been 
treated  classically.  In  classical  theory  spontaneous  and  stimulated  emission  are  two 
different  problems.  Spontaneous  emission  is  obtained  from  Maxwell’s  equations,  treat¬ 
ing  the  electron  beam  as  a  fixed  current  source.  To  obtain  stimulated  emission,  how¬ 
ever,  we  must  consider  the  effect  of  the  electromagnetic  wave  on  the  electron  beam. 


This  can  be  done  using  either  a  collective  approach^  (Boltzmann  equation)  or  a  single 
particle  approach^.  There  is  some  discrepancy  in  the  two  results. 


In  this  paper  we  will  present  a  quantum  mechanical  treatment  of  resonant  transi¬ 
tion  radiation.  In  the  quantum  mechanical  approach,  the  electrons  and  the  electromag¬ 
netic  wave  are  treated  simultaneously  and  both  spontaneous  and  stimulated  emission 
come  out  from  the  same  formalism.  A  simple  fundamental  relation  is  established 
between  the  spontaneous  emission  rate  and  the  stimulated  emission  gain.  The  spon¬ 
taneous  emission  intensity  obtained  quantum  mechanically  agrees  exactly  with  the  clas¬ 
sical  results^.  However,  the  stimulated  emission  gain  obtained  quantum  mechanically 
differs  from  both  Refs.  4  and  5. 

In  the  quantum  mechanical  theory  we  consider  the  interaction  between  the  elec¬ 
tron  beam  and  the  electromagnetic  normal  modes  of  the  periodic  medium.  These  nor¬ 
mal  modes  consist  of  an  infinitu  number  of  spatial  harmonics  with  wavequmbors 

ko  +  (n  =  integer,  kg  =  2ir/\'  =  X/vC  whose  amplitudes  depend  on  the 

Fourier  components  of  E(z).  The  spatial  harmonics  have  phase  velocities  less  than  that 
of  the  fundamental  (n  =  0)  by  the  factor  (1  +  nX /P).  Consequently  it  is  possible  for 
the  electron  beam  to  Cerenkov  radiate  into  the  higher  spatial  harmonics  even  if  v<c. 
Resonant  transition  radiation  can  equivalently  be  viewed  as  a  process  of  Cerenkov 
emission  into  the  higher  spatial  harmonics.  The  results  for  both  spontaneous  and 
stimulated  resonant  transition  radiation  are  the  same  as  those  for  ordinary  Cerenkov 
radiation  multiplied  by  the  squared  amplitude  of  the  n^^  spatial  harmonic  of  an 
appropriate  field  component. 

Section  11  gives  a  brief  summary  of  the  classical  results  for  spontaneous  resonant 
transition  radiation.  In  Section  III  we  describe  the  quantum  theory  of  ordinary  Ceren¬ 
kov  radiation;  both  the  spontaneous  emission  rate  and  the  stimulated  emission  gain  are 
discussed  and  a  simple  fundamental  relationship  is  established  between  them.  These 
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results  are  modified  in  Section  V  for  resonant  transition  radiation  using  the  spatial  har¬ 
monic  amplitudes  obtained  in  Section  IV. 


n.  Classical  Theory: 

The  energy  Uq  radiated  in  a  frequency  range  du;  in  a  solid  angle  dfl  by  an  electron 
traversing  an  interface  between  two  media  with  dielectric  constants  Gi  ^nd  G2  given 


^  ^0  _  e^  sin‘g  A£  ”  ^2  1-// cos&-^~ 

dwdn  4jr^  €  (l-/j’cosfl)(l-/3^^cos^tf) 


where 


-  I  61  €2!  /€o 

€=(€,+  G2/2G0 

Go  -  permittivity  of  vacuum 
13  =  v/c 

Using  Equation  (1),  we  can  rewrite  Equation  (2)  as, 

^^^0  _  e“  sin^g  AG  ^  {nqr^ )  *  <2) 

dudfl  4jrGocvC  4?:^  €  n^q^[2-nq^) 

where  q  =  X/P  vC 

Id  a  periodic  medium  with  multiple  interfaces  the  total  radiated  energy  U  is  given 


•  "  ^  *_•  *  •  ”  ’’m  • 


-■N  .'V.  - 


by 


(3) 


d2u  d-Uo  4  sin-(^(',/2?)sin'(^N/2) 

■  .  •*"  -  -  •  —  -  —  -  ■  -  — 

du/dfl  dcjdn  sin’“(^/2) 

where  is  the  length  of  one  of  the  layers,  2N  is  the  number  of  layers  and 

‘’jrd  1  r=. 

^  =  ■^^(“7  ~  >/tcosfl).  If  N  is  sufficiently  large,  the  radiation  shows  very  narrow  spec- 
X  p 

tral  peaks  around  values  of  0  for  which  ^  is  an  integer  multiple  of  27r.  The  location  of 
the  peaks  is  given  by 

^  =  2t  •  n 


or, 

j  -  yJecosB  =  -5^ 

which  is  identical  to  Equation  (1).  The  narrow  spectral  peaks  can  be  approximated  by 
delta  functions  to  give  for  a  particular  value  of  n. 


d'-U 


d2U„ 


du;dn  du;dn 


SrrNsin^ 


nrrf, 


6(e-2nrr) 


,  nrt,  , 
—  4Nqs,Q  —6 


cos?  +  nq  - 


(4) 


Integrating  over  the  solid  angle  we  get  the  energy  U  radiated  into  the  frequency  inter¬ 
val  do;  for  a  particular  value  of  n. 

12 


dU  . 

ag 

nq-^ 

da; 

^  SIO  V 

4ff6oCv€ 

G 

rr  9 

n^q^(2-nq)9^  j 

(5) 


where  cos?  =  —  nq 

Here  we  have  used  Equation  (3)  for  d^Uo/dwdfl  and  assumed  that  it  varies  slowly  com¬ 
pared  to  the  delta  function.  Equation  (5)  gives  the  total  energy  radiated  by  one  elec¬ 
tron;  multiplying  by  J/e  (J  =  current  density)  we  get  the  radiated  power  per  I  unit 


eJsin"tf  2N 


dcj  4»‘GocVt  G 


n2q2(2  -  )  | 


m.  Quantum  Theory  of  Cerenkov  Radiation: 

In  this  section  we  will  describe  spontaneous  and  stimulated  Cerenkov  radiation 
using  a  quantum  mechanical  formalism.  As  we  have  mentioned  in  the  introduction,  the 
results  for  Cerenkov  radiation  are  readily  adapted  to  resonant  transition  radiation  sim¬ 
ply  by  multiplying  with  the  squared  amplitude  of  the  spatial  harmonic  of  an  appropri¬ 
ate  field  component;  this  is  done  in  Sections  IV  and  V. 

Consider  a  free  non-relativistic  electron  (relativistic  effects  are  incorporated  later 
with  a  simple  modification)  interacting  with  a  radiation  field. 

Ho  =  ^  +  E  (7a) 


Hint  =  EKtX,.  + 


K  -  -  /  P  -  g  L  p-iC*R  + 

V  2m  P" 


where 


k,^'  represent  the  wavevector  and  polarization  of  the  photon  mode, 
Wfc  =  ck/v€ 

V  is  the  volume  of  normalization, 
p  is  the  momentum  operator  for  the  electron, 
m  is  the  electron  mass, 

R  is  the  electronic  position  operator. 
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The  initial  state  1 1>  has  the  photon  modes  in  harmonic  oscillator  states  |  and 

the  electron  in  the  state  |^i>  [  =  (Vr'/^exp  i(]5i*R  -  Ejt/h),  E-,  =  t?;9-^/2mJ. 

ii>  =l^i>  n  (81 

k,i/ 

We  consider  transitions  to  a  final  state  |  F>  with  the  electron  in  state  |^f>  and  one 
more  (or  one  less)  photon  in  mode  k,i/  corresponding  to  photon  emission  (or  absorp¬ 
tion).  The  first-order  matrix  element  M  for  thb  transition  is  given  by 


M  +  1  (emission)  (9a) 

=  \/^u  (absorption)  (9b) 

Assuming  that  the  length  of  the  interaction  region  is  L,,Ly,L,  in  the  x,y,z  directions 
and  that  the  time  of  interaction  is  T,  we  get  from  Equation  (9), 


I  Ml 


2m 


M 


LxLyL,*T 

V 


where 


sin(nxLx/2)  sin(nyLy/2) 

{n,L,/2)  (nyLy/2)  “ 


sin(n,L,/2) 

(n,L,/2) 


=  (^i  ^)x,y.i 


(10) 


n  =  (Ej-Ef  TikJk) 


and  the  upper  (lower)  sign  is  for  emission  (absorption). 

The  total  transition  probability  is  given  by  |  M|  *  integrated  over  the  possible  final  elec¬ 
tronic  states  d^f.  If  L,,Ly,L,  are  sufficiently  large  we  can  treat  the  sin^x/x^  factors  as 
delta  functions  in  this  integration. 


(11) 


±Anp  =  (nj  +  -  ±  -) 


L,LyL,  e¥/3i^sin^g 
2  ~  2'  V  2EeoVm2lkJk 


T2 


sinnT/2  F 
QT/2  ) 


The  electron  couples  only  to  the  photons  which  are  polarized  in  the  plane  containing 
and  Ic  (Fig.  2).  To  account  for  relativistic  effects  we  should  use  the  Dirac  equation 
rather  than  the  Schrodinger  equation  for  the  electrons.  However,  since  electron  spin 
does  not  play  any  role  in  Cerenkov  emission,  we  can  use  the  Klein-Gordon  equation^’*. 
It  is  shown  in  Ref.  7  that  using  the  Klein-Gordon  equation  we  get  m(^^  +  'tf)/2  in  place 
of  m  in  the  expression  for  |  M|  (Equation  10)  where 


E  =  Tfmc* 

'T  = 


(12) 


Consequently  Equation  (11)  is  modified  to 


±  Aiif  =(de+|-±|) 


L^LyL,  e^v^sin^ff 

27i 

V  ^€oVlk^k 

7i  +  7f 

(13) 


where  we  have  used  the  relation  v  = ’h/Jj/m'Yi.  Using  the  momentum  conservation  con¬ 
dition  (n,  y  ;,  =  0)  and  the  relativistic  energy-momentum  relation  (E*  =  p*c^+m*c^)  we 
can  write. 


n  =  Tkvcos® 


J2l_ 

+ 


Twk  + 


m(Tf-,  +  7j) 


=  Tkv 


27i 


7i  +  7f 


1  7f  +  7i  ^  k 
COS0 - 

2Tfi  2/?i 


(14) 


Usually  the  energy  hu/jj  of  the  emitted  photon  is  a  very  small  fraction  of  the  electron 
energy  so  that  Ef  as  Ej  and  -7^  ai  7;  =  7.  Also  we  can  assume  L^LyL,  =  V  so  that  the 
wavefunctions  are  normalized  to  the  interaction  volume.  With  these  assumptions  we 
can  simplify  Equations  (13)  and  (14), 


I  »  _  I  1  .  1  e^v^sin^tf  ,,,2  sin  nT/2 


where  0  =  Tkvcosfl  ±  U|j  +  hkV2Tfm 


Tkv  costf  — ^  T 


(15b) 


The  total  spontaneous  emission  U,  by  an  electron  is  obtained  by  integrating  Equa¬ 
tion  (15a)  (taking  the  upper  sign  with  nj^  =  0)  over  all  photon  modes  H.  If  the  time  of 
interaction  T  is  long  enough  we  can  treat  the  sin^/x^  as  a  delta  function  in  this 
integration  to  get 

^  =  T  —  ~  sinH  (16) 

4jr6oC^ 

To  obtain  the  gain  per  pass  we  need  the  difference  between  stimulated  emission 
and  absorption  which  arises  from  the  small  difference  in  0  for  the  two  cases  (Equation 


_  e“c>;v^sin^g 
"f  4€Q7mc*V 


T2  .p 


where 


^  _  d  sinV  _  16  /  i 

~  TT  “7^  ~  ~T  (maximum) 
d^  I  jr* 


<l>  =  nT/2 

Equations  (16)  and  (17)  give  the  spontaneously  radiated  energy  in  a  frequency  interval 
u  to  M  +  du;  and  the  gain  per  pass  at  a  frequency  u>  due  to  a  single  electron  traveling  at 
a  velocity  v.  For  a  beam  of  electrons  with  a  current  density 
J(=  env  ,  n  =  electron  density)  we  should  multiply  Equation  (16)  by  J/e  to  get  the 
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intensity,  I(W/in^)  and  Equation  (17)  by  nV  to  get  the  total  gain,  F. 


dl  —  /p  .  eJwvsin^l? 

4ir€oC^ 

(18a) 

P  ^  ^  .  eJwvsin^tf  ,  p 

4€o'rnic^ 

(18b) 

Equations  (18a)  and  (18b)  apply  to  ordinary  Cerenkov  radiation;  it  is  shown  in 
Sections  IV  and  V  that  the  results  for  resonant  transition  radiation  are  identical  except 
for  a  multiplying  factor  depending  on  the  relative  amplitude  of  the  appropriate  spatial 
harmonic.  From  Equations  (18a)  and  (18b)  we  can  derive  a  simple  relationship  between 
the  spontaneous  emission  intensity  and  the  stimulated  emission  gain  which  will  hold 
true  even  for  resonant  transition  radiation. 

dl 

r  =  —  ‘(FTr)* — m  (19) 

do;  7 

rV.  Electromagnetic  Waves  in  a  Periodic  Medium: 

Consider  an  electromagnetic  wave  propagating  at  an  angle  9  to  the  z>direction 
which  is  perpendicular  to  the  plane  of  the  layers  (Fig.  2).  The  coupling  of  the  wave  to 
the  electrons  is  proportional  to  the  component  of  the  electron  momentum  along  the 
direction  of  the  electric  field  (Equation  10);  consequently  the  electrons  couple  only  to 
the  waves  polarized  in  the  plane  containing  and  Ic.  For  this  polarization  the  field 
components  are  E,(,E,  and  H^;  the  component  that  determines  the  coupling  strength  is 
E,  since  the  electron  current  is  along  z.  By  Floquet’s  Theorem  E,  can  be  written  as 

E,  =  E,oexp(i(lc-T-wt)]  1  +  ^Pn^xp  i(-^^  + 


(20) 


Starting  from  Maxwell’s  equations,  assuming  an  arbitrary  periodic  dielectric  con¬ 
stant  6(z),  we  can  show  that  for  small  A€^, 


_  a„(nq-^),^ 

Pn  -  — — l2-nq^] 


where 


€(z)  =  €  +  S  *n  +  V'n)- 

n=l 


V.  Quantum  Theory  of  Resonant  Transition  Radiation: 

The  results  (Equations  18)  in  Section  III  have  to  be  multiplied  by  |p„|^  in  order  to 
get  the  results  for  resonant  transition  radiation. 

dl  _^ej(jvsin^g  {^n|^  dq  -  .ggv 

4jr€oC*  **  nq(2  -  nq^ ) 

Assuming  6(z)  has  a  rectangular  form  we  can  readily  calculate  its  Fourier  coefficients 


_  2  Ae  .  "’^1 

■  nx  €  9 


Using  Equation  (23)  in  Equation  (22)  we  get 


dl  _  ^  eJwvsin^fl  A€  .  2  JI.’IIl  n<)  -  . 

dw  4r€QC*  €  ^  n*q*(2  -  nq// ) 


_  eJsin^tf  2Nq  A€  ^  .  9“’^! 

--  —  —  —  ■  ■  *  ■  ■■  —  cin^  ■  ■■ 


4jr6oCV^  ^ 


^  n  V(2  “  nq^ ) 


using  the  relation  VT  =  Nf.  This  is  in  agreement  with  the  classical  result  quoted  ear¬ 
lier  (Eq.  6).  The  stimulated  emission  gain  F  is  readily  obtained  using  the  relationship 

between  F  and  derived  earlier  (Eq.  19)  for  Cerenkov  radiation.  Since  both 
da/ 


quantities  are  multiplied  by  |  p„|  ^  (for  resonant  transition  radiation)  the  relationship  is 
still  valid. 

rV.  ConclttsionB 
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Far>Ultrsviolet  and  X-ray  Emission 
and  Amplifleation  by  Non-relativistic  Electron  Beams 
Traversing  a  Superlattice 

A.  E.  Kaplan  and  S.  Datta 
School  of  Electrical  Engineering 
Purdue  University 
West  Lafayette,  IN  47907 


Abstract 

High-energy  electrons  emit  resonant  electromagnetic  radiation  when  pass¬ 
ing  through  a  spatially  periodic  medium.  It  is  conventionally  assumed  that 
ultra-relativistic  electron  beams  are  required  to  obtain  significant  emission.  We 
demonstrate  theoretically  the  feasibility  of  exploiting  solid-state  superlattices 
with  short  periods  to  obtain  both  spontaneous  and  stimulated  emission  in  the 
far-ultraviolet  and  soft  X-ray  range  using  non-relativistic  beams. 
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Fast-moving  electrons  emit  electromagnetic  waves  when  moving  from  one 
medium  into  another  with  a  different  dielectric  constant  {1-3].  This  is  known 
as  transition  radiation  and  was  predicted  by  Ginzburg  and  Frank  [Ij.  In  a  spa¬ 
tially  periodic  medium  the  waves  emitted  at  different  interfaces  interfere  so 
that  a  resonant  emission  is  obtained  when  the  following  condition  is  satisfied 
(2-4j:  ^ 

V?  costf  =  c/v  -  nX/f  (1) 

where  X  is  the  wavelength  of  radiation,  9  is  the  period  of  the  spatially  varying 
dielectric  constant  e  of  medium  (it  is  conventionally  assumed  that  the  varia¬ 
tions  are  very  small),  v  is  velocity  of  the  electrons  (assumed  normal  to  the 
interfaces),  $  is  the  angle  between  th^  direction  of  wave  propagation  and  elec¬ 
tron  motion,  n  is  an  integer,  and  ?  is  a  “mean”  e.  This  condition  is  readily 
derived  by  requiring  that  the  waves  emitted  at  different  interfaces  interfere 
constructively  at  a  distant  point.  Usually  the  period  0  >>  X,  so  that  ultra- 
relativistic  beams  (v/c  ^  1)  are  required  in  order  to  satisfy  Eq(l)  for  real  0. 
Recently  the  possibility  of  stimulated  resonance  radiation  of  ultra-relativistic 
('^  50  Gev)  electrons  traveling  through  a  stack  of  metal  foik  was  considered  [4j, 
with  J  ~  7  cm. 


In  this  paper,  we  demonstrate  the  feasibility  of  using  non-relativistic  elec¬ 
tron  beams  in  order  to  attain  both  spontaneous  and  stimulated  emission  in  the 
ultraviolet  and  X-ray  range  using  solid-state  superlattices  with  9  ~  100  A  so 
that  9/a\  ~  1.  We  show  that  the  wavelength  of  resonant  radiation  and  the 
required  energy  of  electrons  are  determined  by  the  parameter  Q  =  nXp/J, 
where  Xp  is  a  “mean”  plasma  wavelength  of  the  medium.  If  nXp  <  <  9  (i.e.  Q 
<<  1),  as  is  assumed  in  all  previous  work  [2-4],  then  the  wavelength  of  the 
resonant  radiation  has  an  order  of  magnitude  of  X  X-Q,  and  the  kinetic 
(dimensionless)  energy  of  the  electrons  eU/mc^  must  exceed  the  critical  amount 
~1/\/Q^-^  >>  1  (^  <  Q)  which  constitutes  the  use  of  ultra-relativistic 
beams.  On  the  contrary,  if  the  period  of  the  spatially  periodic  medium  is 
chosen  so  that  Q  >>1  (i.e.  period  9  is  much  shorter  than  the  plasma 
wavelength,  e.g.  9  ~  50*200  A  ),  the  wavelength  of  resonant  radiation  becomes 
of  the  order  of  ~  Xp/Q  =  9/n,  and  the  critical  kinetic  energy  of  the  beam  turns 
out  to  be  extremely  small;  eU/mc^  ~  1/2Q^,  which  may  be  less  than  a  kilovolt 
even  for  very  short  wavelen^h  radiation.  The  advantages  of  the  proposed 
method  are:  (1)  the  frequency  of  radiation  can  be  easily  tuned  in  a  very  wide 
range  by  simply  varying  accelerating  potential  of  beam  (which  is  very  hard  to 
do  with  ultra-relativistic  beams)  (2)  the  range  of  the  possible  angles  of  the  wave 
propagation  k  almost  unlimited,  and  (3)  the  cost  of  equipment  and  energy 
required  for  experiments  with  non-relativistic  beams  k  insignificant  compared 
to  large  accelerating  machines.  This  last  consideration  k  perhaps  the  most 
important. 


The  main  requirements  for  non-relativktic  short-wavelength  resonant  radi¬ 
ation  k  a  periodic  medium  with  a  very  short  spatial  period.  Fortunately,  the 
development  of  molecular  beam  epitaxy  IMBE)  and  other  techniques  in  recent 
years  has  made  it  possible  to  grow  very  tnin  films  (~  100  A  and  less)  with  pre* 
eke  boundaries.  Periodic  structures  composed  of  thin  films  of  different  materi- 
ak,  called  superlattices,  have  also  been  fabricated  (5|.  Using  these  superlat¬ 
tices,  non-relativktic  electron  beams  with  energies  20-200  KeV  can  be  used  to 
generate  radiation  of  wavelength  100-200  A  and  less. 

We  will  first  obtain  j^e  reson^t  wavelength  of  radiation  from  Eq(l)  not¬ 
ing  that  [3]  fv€  =fiv^- 
thicknesses  and  die' 


®i\Ai  hy/h>  where  f|2  and  €|2  are  respectively  the 
ielectric  constants  of  alternating'  layers  forming  the 


superlattice,  (I  =(^i+02  is  its  spatial  period.  For  short  wavelengths 
(X*  <<  Xi*2)i  *1,2  ~  where  Xjo  are  the  plasma  wavelengths  of  the 

two  materi’ab  forming  the  superlattice,  (Xf2  =  4rmcVe^Ni2  where  N*  Is  the 
density  of  electrons).  Thus,  the  mean  dielectric  constant  may  be  written  as 
?  1  -  where  Xp*  =  Substituting  this  into  (1)  and 

solving  it  for  X,  one  gets  the  resonant  wavelengths: 


X  =X„ 


onant  waveiei 
fi  +  costf 

t/3  Q 


where  Q  =  nXp/J;  =  v/c;  Tf  =  (l-/3^)“‘,  and  =  [H-(Q*-sin*0)~']'/^  is  the 
critical  energy  reouired  for  the  excitation  of  resonant  radiation.  For  Q  <<  1, 
X  ^  Xp(Q  ±  y  lcT~'f^)-  Here  Tfc,  -  >>  1.  such  that  only  an 

ultra-relativistic  beam  can  excite  radiation.  On  the  other  hand,  when  Q  >> 
1,  the  critical  kinetic  energy  turns  out  to  be  extremely  low, 
(eU/mc*)„  =  T(cr“l  ^  1/2Q^,  which  is  less  than  500  eV  for  all  conventional 
materials  if  0  ~  100  A.  For  sufficiently  higher  (but  still  non-relativistic)  ener¬ 
gies  eU,  Eq(2)  gives  simply 

X  a  — (-r-cosfl)  2s!  — (\/ mc*/2eU-costf);  eU  <<  me*  =  0.51  MeV.  (3) 
n  p  n 

For  instance,  if  J=100  A,  n  =  1,  eU  =  75  KeV,  and  $  =  45°,  one  has 
X  =  113.6  A;  for  n  =  10,  X  =  11.3  A. 

The  resonant  radiation  (which  may  be  regarded  as  spontaneous  emission  in 
the  system)  can  provide  a  narrowband  source  of  radiation.  A  single  electron 
traversing  multilayer  structure  with  N  layers  radiates  energy  I  in  a  solid  angle 
dn  in  the  fremiency  interval  between  u  and  cj+do;,  given  by  [1-4] 

dH/dudn  =  (d*I„/du/dn)  •  4sin*(e(l,/f)sin*((N/2)/sin2e  (4) 

where  ^  =  (— -\/t’cos0)jrt?/X,  and  d^/dfldw  is  a  radiation  produced  by  a  single 

interface.  According  to  Ginsburg-Frank  theory  [l-3j.  for  non-rel^ivistic  elec¬ 
trons  (^  <<  1)  and  small  variations  of  c  (|  Aej  =  |  C]~e2|  distri¬ 

bution  of  single-interface  radiation  is  given  by 

d  VdHdw  =  e*/9*|Ac(a;)j*sin*tf/4jr*c  (5) 

If  the  number  N  of  layers  is  sufficiently  large,  (N  >>  [c/Acj ),  Eq(4)  provides 
for  very  narrow  spectral  peaks  of  radiation  for  each  particular  angle  $  (with 
central  wavelength  determined  by  Eq(3)),  which  also  implies  that  any  frequency 
is  radiated  in  a  very  narrow  intervals  of  angles.  Noting  that 
Ac  =  X*(Xf*-X2*),  and  integrating  Eq(4)  over  (j  and  (with  dfl  —  2irsin0d9), 
one  gets  the  total  radiation  in  each  order  n 

I  af  16e*Lj*(X,-*-X2-2)*sin*(njrf  ,/J  )/3^u\  (6) 

(where  L  =  Nf/2  is  the  total  thickness  of  the  structure)  with  the  wavelengths 

of  radiation  being  in  the  range  — (“-!)  <  X  <  -(-^  +  1).  The  total  energy 

n  /j  op 

of  radiation  increases  as  speed  decreases.  In  order  to  calculate  the  power  of 
resonant  radiation  emitted  by  an  electron  beam  with  an  electrical  current  J, 
one  must  multiply  Eqs(4)-(6)  by  J/e.  If  ^  =  100  k,9i  =  f2  ”  L  =  1  fim, 
eU  =  75KeV,  J  =  linA,  and  X]  a  400  A  (e.g.  Zn,  Cu,  Ag  or  Au),  X2  a  800  A 
(e.g.  Si  or  Ge,  see  [6j),  the  system  can  provide  a  radiation  of  first  harmonic  (n 
=  1)  with  a  total  power  ~  0.33  mW  and  a  mean  wavelength  ~  200  A. 

We  will  derive  now  an  amplification  caused  by  the  stimulated  emission. 
This  effect  may  be  viewed  in  the  following  way.  An  EM  wave  having  a  wave 


vector  component  k,  =  cos0  along  the  axis  z  (which  coincides  with  the  elec¬ 
tron  trajectory)  produces  the  higher  order  spatial  harmonics  with 
k,  =  k,  ^  2m/f  which  is  due  to  the  periodicity  of  medium.  The  phase  velo- 

0 

cities  ^  these  harmonics  along  the  axis  x  are,  therefore, 
v„  =  c/\/^costf±27m/fko).  If  the  resonant  condition  (1)  for  X  b  fulfilled,  one 
of  these  phase  velocities  coincides  with  the  speed  of  the  electron  that  results  in 
an  exchange  of  energy  between  the  EM  wave  and  the  electron.  For  some  fre¬ 
quencies  in  the  neighborhood  of  resonance,  the  electron  loses  ener^  to  the  EM 
wave;  this  results  in  a  coherent  gain  of  the  wave,  or  stimulated  emission. 

Essentially,  this  resembles  a  common  mechanism  of  amplification  for  many 
kinds  of  microwave  devices  based  on  .the  interaction  of  electrons  with  “slow” 
EM  wave.  The  important  point  is  to  find  the  intensities  of  the  resonant  spatial 
harmonics  of  the  field.  In  all  the  previous  work  on  resonant  radiation  [2-4]  it  is 
assumed  that  X  <  <  (  (  which  is  always  valid  in  the  ultra-relativistic  case,  see 
the  introductory  section).  This  allows  one  to  use  the  WKB  approximation. 
This  approximation  is  not  valid  in  our  case  since  X  may  be  of  order  or  longer 
that  9.  Instead  we  will  find  a  solution  of  the  exact  wave  equation  (with 
periodic  parameters)  based  on  the  assumption  of  smallness  of  variations  of  sus¬ 
ceptibility  (i.e.,  I  Ar/rj  <<1.  which  is  always  true  for  short  wavelengths);  no 
assumption  is  made  regarding  the  ratio  \/9 .  Furthermore,  the  spatial  variation 
of  c(z)  is  usually  approximated  by  a  cosine  function  [2-4j.  In  this  approach,  the 
relative  amplitude  /)„  of  n^  harmonic  of  the  EM  field  is  (Af/t)",  so  that 
for  small  Af  is  negligible  for  all  but  the  smallest  n  (=  1).  In  our  approach, 
we  can  treat  any  arbitrary  function  of  c(z),  in  particular  the  true  rectangular 
function.  We  show  that  falls  off  algebraically  like  the  Fourier  coefficients  of 
c(z).  Significant  radiation  is  expected  even  for  large  n  provided  the  interfaces 
are  sharp  enough.  In  this  letter,  we  approach  the  problem  using  a  single¬ 
particle  picture  which  provides  direct  insight  into  the  mechanism  of  the 
electron-EM  wave  interaction.  The  problem  can  also  be  treated  using  either 
the  Boltzmann  equation  [4]  or  a  quantum  mechanical  formalism;  we  plan  to 
address  these  aspects  in  a  subsequent  publication  [7|. 

We  consider  the  exact  Maxwell  equation  for  the  EM  field  with  e(z)  being 
an  arbitrary  periodic  function  in  z.  We  assume  a  plane  wave;  it  can  be  shown 
that  only  the  EM  wave  with  its  electric  field  E  polarized  in  the  plane  of 
incidence  (i.e.  plane  x,z)  may  be  amplified  by  the  beamJSj.  By  virtue  of  the 
Floquet’s  theorem  for  wave  equations  with  periodic  coefficients  [9],  any  com¬ 
ponent  of  the  EM  field  can_be  written  as  a  sum  of  spatial  harmonics: 

u  =  Uoexp(j^--F-jwt)|l+  X)p„exp(2jnjrz/P  +  j0„)];  (7) 

where  ?*T  =  k^xsintf+koZcos^,  and  is  the  amplitude  of  the  n^  spatial  har¬ 
monics.  We  make  the  conventional  assumption  that  there  is  no  retroreflection, 
which  is  valid  if  N|  Ac/cj^  <<  and  |  Ac/eJ  <<  cos0.  This  assumption  is 
strictly  true  in  the  vicinity  of  0  =  45  *  (see  [8|).  Substituting  the  EM  field  in 
the  (7)  into  the  Maxwell  equations,  collecting  the  terms  with 

exp(jk?+2jn}rz/P)  for  each  particular  n  and  retaining  only  terms  that  are  first 
order  in  a„^<co80),  where  the  a„'s  are  the  Fourie  coefficients  of  €(z): 

e(z)  =  c  +  S  an^os(2nirz/P  +V'„),  one  gets  the  amplitudes  of  the  spatial 
n=l 

harmonics  of  nonvanishing  con^iunts  of  electric  and  magnetic  fields 

(E„Hy,E,):  E,  =  EgCos^;  H„  =  EoVt;  E,  =  -E^sintf,  and 

'  O  '0  o 


an/2 

1  +qsin*d/costf 

q(2costf+q) 

1+qcostf 

1-qcos^q^  , 

^vbere  is  the  ampli^de  of  the  principal  harmonic  of  total  electric  field,  and 
q  =  2jm/Jko  =  Xn/PVc.  Further  calculations  are  based  on  the  conventional 
model  of  energy  exchange  between  the  EM  field  and  an  electron  which  is  used, 
e.g.  in  the  theory  of  free-electron  lasers  (see  e.g.  (10]).  From  the  Lorentz  equa¬ 
tion  me  d(^'7)/dt  =  e(E5  +  l^f!j).  onp  gets  the  equation  for  the  energy 
6  =  Tfmc*  of  electron  _  -  . 

dS/dt  =  ec{M)  =  ec(^E+EA^+ME);  (9) 

where  S  =  S[itt),tj  is  the  field  at  the  instantaneous  location  of  the  electron,  ^ 
and  S  are  unperturbed  v^tors,  is  a  small  perturbation  of  electron  velocity 
due  to  interaction,  and  AE  is  a  small  perturbation  of  the  field  seen  by  the  elec¬ 
tron  due  its  spatial  displacement  in  respect  to  the  unperturbed  trajectory,  i.e. 

AE  =  Az(t);Az  =  c  /  A/9,dt  (10) 

ax  Q 

For  the  assumed  polarization  of  the  field,  it  follows  from  the  Lorentz  equation 
that  ^  ^ 

I  /  E,dt  (11) 

Ifinc  0  -f  me  0 

In  Eqs  (9)‘(11)  one  has  to  take  into  account  only  that  particular  n^  component 
of  the  wave  which  is  “resonant”  to  the  speed  of  electron,  i.e.  that  one  with 
I  4-n;r|  <<  1.  After  substituting  z  =  c^t  and  the  amplitudes  (8)  of  the 
proper  resonant  harmonic  of  E,,Hy,  and  E,,  into  (9)*(11),  integrating  over  the 
temporal  interval  [0,r  =  L//7cL  wnere  r  is  a  time  for  an  electron  to  pass 
through  the  superlattice,  one  has  to  average  the  result  over  all  the  possible 
phases  <t>„  of  the  relevant  field  harmonic.  We  denote  this  operation  by  angle 
brackets.  Note  that  the  term  <^E>  in  (9)  vanishes,  i.e.  the  stimulated  emis¬ 
sion  is  only  due  to  changes  in  the  electron  motion  caused  by  the  field.  Finally, 
one  gets  the  total  averaged  change  of  the  electron  energy  per  pass: 

,  e%*L*  I  ,  c  t'n 

<Ain>  =K  •  - — -  [p  cose{p  cose-^Vep  )  —  X 

m/rc*Tf  («'„f)’  on  n  u; 

( l-cosi/„r) + “f^p^  sin*fl(-2 + 2cosi/„t+  i/„7sini/„r)| ;  ( 12) 

n 

where  =  w{/?v^cos^+q)-l|  =  /?c(njr^^)/j?  is  s  resonant  factor.  For  some  i/„, 
the  change  of  energy  <Ai>  becomes  negative  which  constitutes  the  gain  of 
EM  field,  <A^em^  ~  -<A6>.  In  the  non-relativistic  case,  the  main  contri¬ 
bution  to  the  change  of  energy  is  due  to  z-components  of  A^  and  aE  i.e.,  in 

(12)  —  (/>*,  Pj,Py)  <<  p*  £s<  (a„/2)*(l-/?co8^)*.  Replacing  the  term 

(-2+2co8*/r+wsini/r)/i/*r^  by  its  negative  extremum  -4/n^  {i/t  at  n),  one  gets 
the  maximal  EM-wave  gain  per  electron  per  pass: 

<AIem>  =  a^e^EoL*  sin*^(l/^-cosff)*/jr*mc*/JX  (13) 

In  order  to  obtain  an  amplification  F  per  pass  in  the  system  bombarded  by  an 
electron  beam  with  the  density  of  electric  current  i(A/cm^),  one  has  to  multiply 


-6- 


f 


<Af  em>  ^y  '/®  divide  by  the  energy  flux  of  incident  EM  wave  per  unity 
area  of  the  interface  EoCostf/2R,  where  R  =  3770  is  the  vacuum  impedance. 
One  has  also  to  take  into  account  that  for  rectangular  form  of 
t(z),  a„/2  =  (Af/nir)  sin(njrfi/(!)  with  Ac  =  X^(Xf*-X2^).  Bearing  in  mind  a 
resonant  condition  (-3),  one  finally  gets  the  maximal  EM  wave  amplification  per 
pass: 

r  =  8;iieRL®sin^(7raJ|/J  )sin^fl/mc^Pjr^cosfl,  (14) 

wtiBrc 

/I  =  r’f^(Xr*-X2"^)^(l/^-cos<l)*n-*  =  X'^r*(Xf2-X-2)2(cos«+nX/e). 

If  X  =  140A,  ^  =  lOOA,  U=h  =  50A,  n  =  1,  L  =  l^m,  X,  as  400  A, 
X^  m  800  A,  11  =  45  “ ,  and  i  =  5x10™  A/cm^  (4]  (i.g.  beam  of  2  ftm  diameter 
with  a  current  ~  1.5x10^  A),  one  gets  an  amplification  F  at  5%  per  pass.  The 
required  speed  of  electrons  is  i3  at  0.474  which  corresponds  to  energy  ell  =  69 
KeV.  For  larger  X  the  amplification  increases  drastically.  With  the  mirrors 
situated  outside  the  superlattice  to  form  a  Fabry-Perot  resonator  to  provide 
feedback,  the  system  becomes  a  short*wave  laser  which  may  transform 
significant  portion  of  energy  of  electron  beam  into  coherent  radiation.  It  is 
obvious  that  the  amplifiers  and  lasers  based  on  the  proposed  principle  should 
work  in  the  short  pulse  regime  of  operation,  with  the  duration  of  current  pulse 
being  determined  by  the  heating,  ionization,  diffusion  of  absorbed  electrons, 
etc. 

In  conclusion,  we  have  demonstrated  the  feasibility  of  generating  far- 
ultraviolet  and  soft  X-ray  radiation  by  electron  beams  with  relatively  low, 
non- relativistic  energies,  traversing  the  solid-state  superlattice  composed  of 
very  thin  periodic  layers.  The  use  of  low  energies  is  a  desirable  feature  as  com¬ 
pared  with  ultra-relativistic  beams.  The  proposed  system  can  be  used  as  a  very 
efficient  noncoherent  source  of  narrowband  radiation,  and,  under  special  condi¬ 
tions,  as  an  amplifier  and  laser. 

This  work  is  supported  by  AFOSR  grants. 
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V/e  consklet  the  cyclotron  resonance  in  a  semiconductor  which  has  a  band  with  a  relativistic  energy- momentum  disper¬ 
sion  and  show  that,  when  subject  to  quasi-resonant  radiation,  the  energy  of  a  free  carrier  in  this  band  and  the  relative  po¬ 
tential  drop  between  the  two  surfaces  of  the  semiconductor  facing  the  radiation  exhibit  bistability  under  classical  single 
electron  apptoximatioiL 


In  a  previous  work  by  one  of  the  authors  [  1  j ,  it 
was  shown  that  cyclotron  resonance  of  relativistic 
electrons  exhibits  bistability  due  to  energy  dependence 
of  the  relativistic  cyclotron  frequency.  In  the  present 
letter,  we  consider  the  possibility  of  bistable  cyclotron 
resonance  in  a  semiconductor  which  has  a  nonparabolic 
band  with  a  relativistic  momentum  dispersion  in  at 
least  part  of  the  first  Brillouin  zone  (2—4] .  By  means 
of  a  classical  analysis,  we  demonstrate  that  the  cyclo¬ 
tron  of  free  carriers  in  the  semiconductor  exhibits 
bistability  if  the  semiconductor  is  immersed  in  a 
homogeneous  magnetic  field  Rq  and  irradiated  by  a 
circularly  polarized  light  propagating  along  Rq, 

The  present  analysis  is  similar  to  that  of  ref.  [  1  j . 
However,  there  are  a  number  of  novel  features  in  the 
present  problem.  First,  the  nonlinearity  of  free  carriers 
in  semiconductors  are  several  orders  of  magnitude  larger 
than  the  nonlinearity  of  relativistic  electrons  in  vac¬ 
uum.  This  lowers  the  threshold  intensities  for  the 
onset  of  hysteretic  jumps  in  semiconductors,  despite 
rapid  thermal  relaxation  which  broadens  the  cyclotron 
resonance.  In  contrast  to  ref.  (1],  the  strength  of  the 
nonlinearity  does  not  allow  one  to  neglect  the  veloc¬ 
ity  terms  higher  than  the  fourth  order  and  the  har¬ 
monics  higher  than  the  first.  We  give  an  exact  result 
in  the  steady  state.  Second,  effective  masses  can  be 


quite  small  in  semiconductors.  For  instance,  in  narrow 
gap  semiconductors  they  are  nearly  two  orders  of 
magnitude  smaller  than  the  bare  electron  mass  in  the 
vicinity  of  the  band  edge  [5].  This  increases  the  cy¬ 
clotron  frequency,  making  it  possible  to  pump  carriers 
with  optical  or  near  optical  frequencies.  For  example, 
if  the  semiconductor  is  n-type  InSb  and  ~  1^0  kC, 
one  can  use  C02-laser  at  10.6  ji.  Third,  the  nonlinear 
cyclotron  resonance  in  semiconductors  is  accompanied 
by  the  appearance  of  a  voltage  drop  between  the  two 
surfaces  of  the  sample  facing  the  radiation.  The  voltage 
drop  arises  from  the  redistribution  of  rotating  carriers 
under  the  influence  of  the  radiation  pressure  and  ex¬ 
hibits  bistability.  Therefore,  the  proposed  effect  could 
be  the  first  known  all-optical  nonlinear  phenomenon 
which  yields  an  opto-electronic  bistability. 

Consider  a  semiconductor  with  a  conduction  band  * ' 
which  has  a  relativistic  momentum  dispersion  as  in 
Kane's  isotropic  two-band  model  [2], 

The  argument  is  easily  generalized  to  holes  of  a  valence 
band  with  a  relativistic  momentum  dispersion. 

This  type  of  semiconductors  have  been  discussed  exten¬ 
sively  for  the  purpose  of  obtaining  a  cyclotron  maser.  See 
rels.  (3,  4). 


0.031-9163/83/0000-0000/S  03.00©  1983  North-HoUand 


30.S 


Volume  9T-A,  number  7 


PHYSICS  LETTERS 


12  September  198j 


Fig.  1.  (p/po)^  versu»///o  for  cjJ  =  lOr:  (a)  ui  =  4r; (b) cu 
=  sr:(c)u)  =  6r;  (d)w  =  7r. 


+  1]  [(6j*/w)2/3  -  1]-2 

^(w/n2<(a;‘/r)2.  (6) 

Of  course,  to  observe  the  cyclotron  resonance,  one 
must  have  F  ^  co* ,  co. 

It  is  seen  from  fig.  1  that,  as  /  is  increased  from 
zero,  the  system  moves  along  the  lower  branch  of  a 
curve  corresponding  to  co  until  /  reaches /2.  At  / 

=  I2,  the  system  makes  a  hysteretic  jump  to  the  stable 


f  ig.  2.  g/go  versus  1/Iq  for  wg  =  lOr.  a,  b,  c  and  d  are  the 
ame  as  in  fig.  I . 


upper  branch  and  moves  along  the  upper  branch  with 
increasing  I.  Prior  to  the  jump  at  I2,  the  momentum  is 
approximately  given  by 

Pi  . 

If  /  is  decreased  now,  the  system  moves  down  the 
upper  branch  past  I2  to  /j .  At  /  =  /j ,  the  system 
jumps  to  the  lower  stable  branch.  Prior  to  the  jump, 
the  momentum  is 

Pi  -  1]*^^  . 

After  the  jump,  the  momentum  is 

Note  that  the  hysteretic  jumps  at  /j  and  1 2  are  con¬ 
tingent  upon  the  charge  distribution  being  in  steady 
state,  at  least  relative  to  cyclotron  transitions.  In 
other  words,  the  charge  relaxation  must  be  fast  such 
that  ToF  <  1 .  Furthermore,  variations  in  /  must  be 
slow  such  that  <1.  Using  the  expression  for 

mobility  n  =  e/m* F  and  the  conductivity  a  =  eNu, 
where  N  is  the  free  carrier  density,  one  can  write 

TqF  =  eom*F2(47re2AO~  ‘  =  (eo/e-)(F/Up)2  <  1. 

Wp  is  the  plasma  frequency.  This  condition  puts  a 
lower  limit  on  the  carrier  density  such  that  *  ’ 

A^>yVj  =eom*F2(47re2)-l  . 

If  yv <  ,  which  corresponds  to  TqF  ^  1 ,  then 

charge  oscillations  may  be  excited  and  the  system 
may  not  reach  a  stable  steady  state. 

There  is  also  an  upper  limit  on  the  free  carrier  den¬ 
sity  N.  This  comes  about  from  the  fact  that  at  low 
temperatures,  band  states  are  filled  up  to  a  Fermi  mo¬ 
mentum  pp  and  cyclotron  transitions  should  excite 
the  electron  to  states  for  which  p  >pp.  One  can  see 
from  the  curves  in  fig.  1  that  the  transition  from  the 
higher  branch  to  the  lower  cannot  occur  if  P3  <pp. 
Therefore,  in  order  to  observe  the  bistable  behavior 
in  full,  one  must  have P3  >pp,  and  since pp 
=  ft(37r2yV)U3  as  temperature  goes  to  zero, 

eom*F(4rre^)~  *  =N^<N 

<N2  =  (37r2h3)-  lp3[(^j*/co)2/3  -1)3/2.  (7) 

m*  here  is  the  appropriate  effective  mass  near  the  Fermi 
surface.  Typically  m*  ~  Am^. 


